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ABSTRACT 


Relativistic  treatments  of  irreversible  processes  in  a  simple 
fluid  have  been  given  by  C.  Eckart  and  by  Landau  and  Lifshitz.  Superficially, 
the  two  formulations  bear  little  resemblance  to  each  other.  The  basis  of 
the  two  theories  has  been  explored  by  kinetic-theoretical  methods  and  it  is 
understood  that  the  difference  between  them  is  due  simply  to  a  different 
choice  of  rest  frame. 

In  the  present  thesis,  we  have  attempted  the  construction  of  a 
general  phenomenological  theory  which  is  not  tied  at  the  outset  to  a 
particular  choice  of  rest  frame,  and  which  will  include  the  Landau-Lif shitz 
and  Eckart  formalisms  as  special  cases.  A  general  formulation  of  this  type 
leads  to  a  number  of  aesthetically  satisfying  results,  for  instance  the 
amalgamation  of  Fick's  law  of  particle  diffusion  and  Fourier's  law  of  heat 
conduction  into  a  single  law,  in  conformity  with  the  relativistic  equivalence 
of  mass  and  energy.  The  extension  to  fluid  mixtures,  and  the  incorporation 
of  electromagnetic  effects  (of  interest  in  plasma  theory)  are  particularly 
simple . 


The  complete  flexibility  in  the  choice  of  rest  frame  also  brings 
with  it  certain  practical  advantages,  since  a  judicious  choice  will  often 
simplify  the  mathematics.  We  have  illustrated  this  (a)  by  giving  new  and 
simple  proofs  of  the  theorems  of  Tolman  and  Klein  on  the  variation  of 
temperature  and  thermal  potential  in  a  fluid  in  thermal  equilibrium  in  a 
static  gravitational  field;  (b)  by  integrating  the  equations  for  one¬ 
dimensional  steady  flow  in  the  presence  of  viscosity  and  heat  conduction 
to  obtain  an  estimate  for  the  thickness  of  a  relativistic  plane  shock  layer. 
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CHAPTER  I 

INTRODUCTION. 

The  body  of  the  thesis  will  be  concerned  with  construction  of  a 
relativistically  invariant  theory  of  the  thermodynamics  of  irreversible 
processes.  The  classical  theory  of  irreversible  processes  (and  also  the 
two  relativistic  theories  of  Eckart  [1]  and  Landau-Lifshitz  [3])  consists 
of  a  linear  model  of  the  irreversible  processes  based  on  the  theory  of 
reversible  processes.  Such  a  theory,  of  course,  can  be  expected  to  apply 
only  to  processes  near  to  equilibrium. 

Before  we  embark  on  the  construction  of  our  theory  (also  a  linear 
approximation),  it  may  be  helpful  to  review  the  classical  theory  of  reversible 
and  irreversible  processes.  While  so  doing  we  shall  pinpoint  the  assumptions 
which  brand  it  as  a  linear  approximation. 

Section  1:  The  Theory  of  Equilibrium  (Reversible)  Thermodynamics. 

In  the  state  of  equilibrium  a  simple  fluid  may  be  described  by 
three  parameters: 


(i) 

u  = 

energy  per  gram 

1.1.1 

(ii) 

V  = 

volume  per  gram 

(iii) 

s  = 

entropy  per  gram 

The  entropy,  may  be  .expressed  in  terms  of  u  and  v  by  an  equation 

of  state 


1.1.2 


s 


s(u,v) 


We  may  define  T  and  P  by 


o 


1.1.3 


Tds  =  du  +  Pdv 


thus  making 


1.1.4 


T  = 


-1 


P  =  (I*)  T  . 

ovu 


If  we  let  dq  (not  a  perfect  differential)  be  the  amount  of  energy 
transferred  from  the  surroundings  to  one  gram  of  the  fluid,  we  may  postulate 


1.1.5 


Tds  =  dq 


It  is  easy  to  see  that  P  is  the  mechanical  pressure,  since 


1.1.6 


dq  =  increase  in  energy  +  work  done 

of  1  gm  of  fluid  by  1  gm  of  fluid 


Section  2:  The  Classical  Theory  of  Non-equilibrium  (irreversible)  Processes. 

We  begin  by  assuming  (Assumption  l)  that  an  equation  of  state 

1.2.1  s  =  s (u, v) 

still  exists  during  the  irreversible  process.  Next  we  assume  that  s  is 
differentiable  (Assumption  2)  so  that  we  may  write 

1.2.2  Tds  =  du  +  Pdv 

If  we  put 

q  =  energy  flux  vector 


1.2.3 


t,.  =  mechanical  stress  tensor 


-  3  - 


we  may  express  the  conservation  of  energy  as 


1.2.4 


du 


dv. 


—  +  vdivq+vt..  3 —  = 

dt  n  ij  ox. 


=  0  . 


•  • 

We  may  define  the  viscous  stress  tensor  t  by 


1.2.5 


i  j  i  j  _~ 

T  =  t  -  P& 


ij  * 


The  equation  of  continuity 


1.2.6 


in 

dt 


=  -  p  div  v 


0  -  v 


enables  us  to  write  1.2.4  as 


1.2.7 


ds  j.  /q \ 

o  dT  +  dlv  (*> 


'1 

<  ^  - 


-T.  .  dv.  . 

LL  ^  +  q  •  V  (|)  . 


t  3x7 

j 


Observe  that  for  an  isolated  volume  v 


1.2.8 


3F  pAdv  = 


V 


dA 

dt 


p  dv  + 


A  (pdv) 


v 


dA 

it  pdv  • 


v 


The  second  law  of  thermodynamics  asserts  that  the  total  entropy 


of  an  isolated  system  is  non-decreasing.  Thus 


-  4  - 


1.2.9 


d_ 

dt 


ps  dv  >  0 


In  view  of  1.2.8,  1.2.7  snd  the  fact  that  q  =  0  on  the  boundary 
of  an  isolated  region,  we  may  put 


1.2.10 


d_ 

dt 


ps  dv  = 


dv. 

T.  .  v— ^  +  q 

ij dx7  4 


.7- 

T 


dv  >  0 


v 


To  put  the  entropy  production  on  a  local  basis  we  will  suppose 
that  entropy  can  "flow".  Thus 


1.2.11 


Entropy  Production  =  Entropy  Increase  +  Entropy  Export 


s  '  n  dA 


=  ^  ///  ps  dv  + 


v 


V 


where  s  is  the  vector  representing  the  "flow"  of  entropy.  The  rate  of 
entropy  production  (cr)  per  unit  volume  is  then  by  1.2.11 


1.2.12 


ds  , .  — * 

a  =  p  —  +  s 


Next,  to  identify  the  rate  of  entropy  production  at  a  point,  we 
assume  (Assumption  3)  that 


q 

T 


1.2.13 


s 


5 


This  identifies  the  local  entropy  production,  for  according  to 
equations  1.2.7,  1.2.13 


1.2.14 


dv 


1  ij  i  ,lx 

a  =  "  T  T  5^  +  *&(?)  * 


If  we  now  assume  that  cr  >  0  (Assumption  4)  ;  that  t1^  and  q* 

^Vi  1 

are  linear  functions  of  and  ^/  (— )  (Assumption  5)  we  may  deduce 

Fourier's  law 


1.2.15 


"+  / 1  \ 

q  =  k.v(-) 


K  >  0 


and  the  form  of  the  viscous  stress  tensor 


1.2.16 


t'J  -  -a 


dv.  dv 


7  *  >  ' 


6v 


b  t— ^  8 

ij 


a,  b  >  0  , 


Section  3:  Discussion  of  the  Assumptions  of  Classical  Irreversible 

Thermodynamics , 

Let  us  set  down  the  assumptions  made  in  Section  2. 

1.  An  equation  of  state  s  «*  s(u,v)  exists  under  non-equilibrium 
conditions . 

2.  The  relation  Tds  =  du  +  Pdv  holds  under  non-equilibrium 


conditions . 


o 


3.  A  linear  relation  exists  between  s  and  q 

4.  A  linear  relation  holds  between  the  "fluxes” 


dv. 


th 


e"  forces"  ^7  Gf);  * 


5.  The  rate  of  entropy  production  is  positive 


q 

T  * 


and 


The  assumptions  of  linearity  in  3?  4,  5  suggest  that  the  theory 
is  only  applicable  near  equilibrium.  If  we  appeal  to  kinetic  theory  (see 
Prigogine  [11])  we  find  that,  for  gases,  only  (5)  can  be  justified  for 
arbitrary  deviations  from  equilibrium  (it  is  I oltzmanns  H  Theorem).  The 
assumptions  1  -  4  can  only  be  justified  for  near-equilibrium  conditions. 
It  appears  then,  that  the  method  of  section  2  supplies  only  a  linear 
model  for  what  is  most  likely  a  non-linear  physical  situation. 


It  is  worth  noting  that  assumption  3  maY  be  replaced  by  the 
weaker  assumption 


1.3.1 


div  s  =  div 


In  any  case,  the  energy  flux  q  is  only  determined  to  within  curl 
A*  }  when  A*  is  any  vector  whose  curl  vanishes  on  the  boundary  of  an  isolated 
system.  After  all,  we  can  only  ascribe  unique  meaning  to  heat  flow  through 
a  closed  surface,  i.e.  to 


n  dA 


fS  closed) 


1.3.2 


S 


-  7  - 


If 

— >  — ►  — > 

1.3*3  q^.  =  q  +  curl  A 


we  have 


1.3.^ 


(S  closed) 


Hence,  from  the  point  of  view  of  measurement f  is  just  as  good 

a  candidate  as  q  for  use  in  the  equation  of  conservation  of  energy.  If 
we  do  use  it,  we  get  Fourier’s  law  in  the  form 

1.3.5  q**  -  K-V(f) 


which  is  different  from  the  result 

q  -  K  V  (|)  • 


However,  Fourier’s  law  for  a  closed  surface  S  is 


Since 


the  result  for  closed  surfaces  (the  only  one  with  unique  physical  meaning) 

— > 

is  the  same  for  either  choice  of  q 


-  8  - 


Section  4:  Relativistic  Formulation. 


The  points  made  in  section  3  do  not  preclude  construction  of 
a  regorous,  self  contained  mathematical  model  of  irreversible  thermodynamics 
based  on  assumptions  1-5.  Such  a  model  of  course,  would  have  its  range 
of  applicability  restricted  to  near  equilibrium  conditions.  It  would  be 
analogous  to  Newtonian  Mechanics  -  an  elegant  model  of  the  real  world  whose 
applicability  is  restricted  to  ’'classical"  conditions.  In  fact,  such 
models  have  been  constructed,  using  the  frame  work  of  relativity,  by  Eckart 
[1]  in  1940  and  later  by  Landau  and  Lifshitz  [3]. 


The  two  theories  suffer  from  a  basic  lack  of  invariance  with 
respect  to  the  rest  frame  of  the  observer.  Eckart  describes  the  behaviour 
of  the  fluid  as  seen  by  an  observer  who  sees  no  net  particle  flux.  Landau- 
Lifshitz  developed  their  theory  from  the  point  of  view  of  an  observer  who 
sees  no  net  energy  flux.  If  we  think  of  Eckart  and  Landau  as  equipped 
with  small  flags,  standing  in  their  respective  reference  frames,  the 
situation  is^ 


* 


Particles 

stationary 


•  & 


Particles  appear 
to  drift 


'Mr 


Flag  sags 

(no  "energy  wind") 


ECKART 


LANDAU 


-  9  - 


Eckart 's  flag  will  flap  in  the  energy  "wind"  but  he  will  see 
the  particles  as  stationary  on  the  average;  Landau's  flag  will  be  limp 
on  its  staff  but  he  will  see  a  net  drift  of  particles  past  him. 

The  results  obtained  by  the  two  men  differ  in  form,  but,  as  we 
shall  see,  are  consistent.  Landau  describes  the  observed  particle  flux  as 
proportional  to  (  -  )  the  gradient  of  thermodynamic  potential;  Eckart 
describes  the  observed  energy  flow  as  proportional  to  (  +  )  the  gradient  of 
thermodynamic  potential.  The  form  of  the  mechanical  stress  tensor  deduced  by 
the  two  is  identical.  Since  the  velocities  on  which  the  stress  tensor  depends 
are  different  for  the  two  observers,  the  result  is  in  fact  at  variance.  We 
will  show  that  the  difference  is  negligible. 

In  the  process  of  demonstrating  the  equivalence  of  the  two  theories, 
we  have  constructed  a  third  theory  which  is  not  tied  at  the  outset  to  the 
rest  frame  of  any  particular  observer.  By  suitable  specialization,  we  may 
obtain  Eckart' s  theory  or  Landau's  theory  as  a  special  case. 


998  XXiw  ori  3od  "bnlw”  Y,Txsn3  sdl  ai  qftX3t  XXlw  3&X3  e  l?sdd8  g 
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CHAPTER  II 

A  GENERALIZED  FORMALISM  FOR  THE 
IRREVERSIBLE  THERMODYNAMICS  OF  A  SIMPLE  FLUID 


Section  1. 


In  this  chapter  we  will  consider  the  thermodynamics  of  a  simple 
fluid  -  one  composed  of  particles  of  only  one  kind.  It  will  be  described 
microscopically  by  a  symmetric  tensor  T^V  and  two  vectors  MU  and  S^  , 
which  we  shall  take  as  given  a  priori.  We  will  from  the  outset  assume  the 
conservation  equations 

2.1.1  T^V  =  0 

2.1.2  M^1  =  0 

and  the  principle  of  increase  of  entropy 

2.1.3  S  ^  >  0  . 

We  will  use  the  framework  of  general  relativity,  and  utilize  a  metric  with 
signature  +  2.  As  for  notation,  a  comma  stands  for  covariant  differentiation, 
Greek  indices  run  from  1  to  4,  and  Latin  indices  run  from  1  to  3*  Although 
it  is  irrelevant  to  a  formal  phenomenological  description  of  the  fluid,  it 
may  be  helpful  to  explain  the  meaning  of  these  quantities  from  a  microscopic 
point  of  view.  To  do  this  consider  an  observer  0  ,  moving  with  4-velocity 
Ua  ,  and  carrying  a  small  imaginary  box  of  dimensions  dx  :  dy  :  dz  whose 
faces  are  perpendicular  to  the  x  :  y  :  z  axes  of  his  local  (Cartesian) 
coordinate  system.  This  observer  will  denote  quantities  measured  in  his  local 
coordinate  system  by  a  star  (*)  ;  in  particular  =  (0, 0,0,1)  .  For 


11 


simplicity  we  will  use  units  wherein  c  =  1  . 

A.  m!1  has  the  property  that  |M^N^|dE  is  the  number  of  world  lines  of 
fluid  particles  which  intersect  a  3“dimensional  subspace  with  volume  dZ  , 
multiplied  by  the  rest  mass  of  one  particle  (N^1  is  the  normal  to  dZ)  . 

The  meaning  of  this  statement  for  0  is  easy  to  see.  Taking  the  subspace 
dZ  to  be  his  box,  *  (0, 0,0,1)  ,  dZ  =  dV  =  dx  dy  dz  and 
M^N^IdZ  =  M4dV  .  Since  his  box  is  part  of  the  hypersurface  t  =  constant, 
M^dV  is  the  rest  mass  of  the  particles  in  his  box  at  time  t  .  On  the  other 
hand,  making  the  subspace  dZ  to  be  the  history  of  the  face  of  the  box 
perpendicular  to  the  x  axis  during  time  dt  ,  we  have 

t  (1,0, 0,0)  ,  dZ  =  dy  dz  dt  , 

and 

|M  lf|dZ  *  M1  dA  dt  . 

Since  our  subspace  is  now  a  portion  of  the  hypersurface  x  =  constant, 
M1  is  the  total  rest  mass  of  the  particles  crossing  unit  area  perpendicular 
to  the  x  axis  in  unit  time.  Similar  results  hold  for  M2  and  .  To 

summarize:  M*  ,  £  ,  M5  are  the  total  rest  masses  of  the  particles  crossing 

4 

unit  areas  perpendicular  to  the  x  axis,  y  axis,  z  axis  in  unit  timej  M 
is  the  density  of  rest  mass. 

u, 

The  significance  of  the  conservation  equation  M  ^  =  0  is  now 
easy  to  see.  Considering  a  small  region  R  at  rest  in  0's  local  coordinate 
system  and  the  fact  that  in  local  cartesian  systems  covariant  differentiation 
reduces  to  ordinary  partial  differentiation, 
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2.1.4 


dV 


Using  Gauss’  theorem  on  the  R.H.S. 

2-!-5  ^  fff  M4  dV  i  -  JJ  M1  n.  dA  . 

R  R 

2.1.5  states  that  the  number  of  particles  is  conserved,  i.e.  the  only  change 
in  the  number  of  particles  inside  the  volume  is  due  to  a  flux  of  particles 
across  the  surface  of  the  region. 


LtV 

B.  1  has  the  following  properties* 

uv  ■)£  1±1± 

1 .  T  U  U  dV  =  T  dV  is  the  total  energy  inside  a  region  of 
I-  * 

volume  dV  . 


2.  5*  +  u'u  )  U  dV  =  T^dV  is  the  i^  component  of  the 

|J.  p.  v 

total  momentum  inside  a  region  of  volume  dV  . 


3. 

of 

xj 

(dA 


T^V(bl  +  U*U  U 

p  p  v  v 

the  force  exerted  on  a  face  of  the  box  perpendicular  to  the 
axis  by  the  material  corresponding  to  smaller  values  of  x^  . 
is  the  area  of  the  face.) 


\  *  i  i 

)dA  =  T  JdA 


is  the  i^  component 


LLV 

The  interpretation  of  the  conservation  equation  T  ^  =0  is  now 
easy  to  see.  Consider  the  total  force  on  the  material  in  the  observers  box, 
The  forces  tending  to  push  it  in  the  direction  of  positive  x1 


T.-  dx.  dx  +  T.0  dx  dx  +  T  dx  dx0 

ll  2  3  i2  1  3  i3  1  ^ 


are 


■ 


t 


-  13  - 


The  forces  tending  to  push  it  in  the  direction  of  negative  x1 


are 


^Til  dT 

(Til  +  — r  dxl)  dx2  dx3  -  (Ti2  +  -jr  dx2>  dx2  dx3 


dT 


(Ti^  +  dx^)  dxx  dx2 


Thus  the  total  force  pushing  the  material  in  the  box  in  the  direction  of 


positive  x  is 


ciT  , 


u 


dV  . 


dT.  „ 

By  virtue  of  =  0 


6T 


2.1.6 


a  t 


14  dV  S 


6t  . 

— ^  dV 


th 


This,  since  T^  dV  is  the  i  component  of  the  total  momentum  in  the  box, 
just  means  that  the  rate  of  change  of  momentum  is  equal  to  the  force  applied. 


ma 


Consider  next  the  flow  of  energy.  In  relativity,  energy  p  has 
ss  ^  ,  an  energy  flow  (j.  has  momentum  ^  ,  and,  conversely, 


momentum  m  implies  an  energy  flow  m  c  .  Consider  now  a  small  region  R 

fixed  with  respect  to  the  observer  0  .  By  virtue  of  T  ^  =  0 

1 11 


dT 


R 


6x 


V- 


dV  =  0 


or,  using  Gauss'  theorem 


-  14  - 


2.1.7 


The  L.H.S.  is  the  time  rate  of  change  of  energy,  the  R.H.S.  (with  c  =  l) 
is  the  influx  of  energy  across  the  boundary. 


of  the 


To  summarise,  the  conservation 
conservation  of  energy  (T 
0). 


MV 

equations  T  ^  *  0  are  expressions 
and  the  conservation  of  momentum 


C.  is  the  entropy  flux  vector.  Since  we  are  considering  a  fluid  with 

irreversible  changes  taking  place  within  it,  entropy  is  being  created  within 
it.  Thus,  the  entropy  outflux  from  a  region,  plus  the  entropy  increase 

within  the  region,  must  be  positive.  Taking  S1  as  the  entropy  flux  vector 

4 

and  S  as  the  entropy  density 


2.1.8 


dA 


+ 


d_ 

dt 


dV  >  0 


or  using  Gauss  theorem  for  the  arbitrary  (small)  region  R 

2.1.9  S  ^  >  0  . 

t  M- 

LLV 

D.  It  is  worth  noting  that  the  timelike  eigenvector  of  T  has  great 
physical  significance.  Define  by 

2.1.10  T^V  E,  =  -  <P  J  E  =  -  1  . 

V 


In  a  coordinate  system  where 

2.1.11  *  [0,0,0, 1]  , 


-  15  - 


we  have 


_(i4  *  LI 

T  =  -  cp  E 


Thus  E^  is  a  unit  vector  parallel  to  the  momentum. 
Section  2;  Reversible  Thermodynamics. 


Every  fluid  (which  is  not  externally  constrained  to  do  otherwise) 
tends  in  time  to  achieve  a  state  of  equilibrium.  This  state  is  characterised 
by  3  properties; 


1.  There  is  no  heat  flow,  and  no  viscous  forces  act.  Taking  0  to  be  an 
observer  who  sees  the  fluid  as  possessing  no  net  momentum,  this  means 
that 

;  E11  t  (0,0,0, 1) 

Note  that  in  2.2.1  we  have  assumed  the  fluid  to  be  isotropic  so 
that  T11  =  T22  =.  T^  =.  P  .  From  2.2.1  it  is  clear  that 

2.2.2  T^V  =  (\i  +  P)EM'  EV  +  Pg^V  . 

This  we  will  take  as  our  first  characteristic  of  a  fluid  in  equilibrium. 

2.  Since  the^e  is  no  heat  flow,  the  energy  flow  vector  is  parallel  to  the 
particle  flow  vector: 


2.2.1 


,pv 


P  0  0  0 
0  P  0  0 
0  0  P  0 

0  0  0  \i 


2.2.3 
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This  we  will  take  as  our  second  characteristic  of  a  fluid  in 
equilibrium. 


3.  There  exists  a 

2.2.4 


functional  relation  between  the  three  variable 
p  =  f(n,p)  . 


(1>  P  j  P  5 


It  follow  from  the  usual  postulates  of  equilibrium  thermodynamics 
that  the  Pfaffian  form 


d (j”)  +  P(n,p)  d(^) 

P  P 

has  an  integrating  factor  --t-1-  ■  -v  which  is  strictly  positive: 

*  v  P  / 

T(n,p)  d  S(n,p)  =  d(£)  +  P  d (— ) 

P  P 

2.2.5 

T(n,p)  >  0 


S  is  defined  by  2.2.5.  Let  us  assume  that  an  observer  who  sees 
no  net  particle  flow  also  sees  no  net  entropy  flow,  i.e.  assume 

2.2.6  =  S  . 


Then  equations  2.2.2,  2,2.3  and  2.2.6,  plus  the  conservation 
equations  2.1.1  and  2.1.2  are  sufficient  to  force  the  conclusion  that  S 
is  constant  for  a  fluid  in  equilibrium: 

From  2.2.2  and  2.1.1 


0  =  E  T^V 

H  >v 


V  /  \  v 

M  „  E  -  (|i  +  P)  E  . 

t  V  f  V 


2.2.7 


»  » 
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From  2.2.3  and  2.1.2 


2.2.8 


Substituting  2.2.8  into  2.2.7  gives 


2.2.9 


[|i 


>v 


E 


=  0 


Recalling  2.2.6,  2.2.5,  2.2.3  and  2.2.2 


2.2.10 


=  S  p  E 


=  T'1  [h 


(~)  P  1  EV  . 
P  >  V 


It  follows  from  2.2.9  that 


2.2.11 


S  =  0 


and  hence  by  2.2.10 


2.2.12 


S  E  =  0  . 


For  a  fluid  in  equilibrium  and  E|a  are  parallel,  so  E^  is  th< 
only  natural  choice  for  the  4-velocity  of  the  fluid.  Thus 


2.2.13 


S  E 


dS 

dT 


=  0 


and  the  change  in  entropy  of  the  fluid  with  respect  to  proper  time  is  zero. 
Tims  we  have  the  important  results  If  the  entropy  of  a  fluid  in  equilibrium 
is  spatially  constant  at  some  initial  time,  then  it  is  universally  constant. 


•  ’ 
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Section  3»  Irreversible  Thermodynamics, 

Here  we  will  consider  a  fluid  far  from  equilibrium.  Since  heat  flow 
will,  in  general,  be  present,  the  vectors  and  will  no  longer  be 

parallel.  There  is, then,  no  natural  choice  for  the  fluid  4-velocity.  We 
will  retain  the  postulates  of  conservation  of  particles,  conservation  of 
energy,  and  conservation  of  momentum  in  the  forms 


2.3.1 

T^V  =  0  ,  M  ^  =  0 

As  dissipative  processes  are  taking  place,  we  shall  have  in  place  of  S^  =  0, 

}  M- 


2.3.2 

SM  >  0 

- 

Let  us  denote  the  limit  of  any  quantity  Q  ,  as  t  approaches  infinity,  by 
.  If  is  again  defined  by 


2.3.3 

T^V  E  =  -  9  E4  ;  E^  E  =  -1  , 

v  p 

it  is  reasonable  to  suppose  that 


2.3.4 

lim  T^V  =  (p  +  P  )EM'  EV  +  P  g^V 

'  *00  00  '  00  00  00 

t— >  00 

2.3.5 

lim  M^*  =  P„ 

t— >  00 

2.3.6 

Pco  -  "  J  • 

Let  us 

now  choose  some  network  of  observers!  i*e*;  define  a 

4-velocity  field 

2.3.7 

U  =  -1  ;  lim  =  E^ 

H  t  -><*> 
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The  second  requirement  assures  us  that  the  network  of  observers 
will  be  those  of  section  2.2  when  equilibrium  has  been  reached. 

Now  let  us  define  the  quantities  p,  S^V  by 


2.3.8 


2.3.9 


2.3.10 


pU^UV  +  {;^uv  +  +  s^v 

It 


0 

0 


The  physical  interpretation  of  p,  ,  SUV  is  easy  to  get  by  considering 
the  equations  8-10  in  the  frame  where  =  (0, 0,0,1)  .  It  is  easy  to 

see  that 


2.3.H 


< 


,n4 


* 


T 


iU 


T 


ij 


Hence 

p  =  energy  density 

_  energy  flux  (or  momentum  ) 
SMV  =  stress  tensor 


measured  by  the  chosen  network  of  observers. 

Invariant  expressions  for  p  ,  ^  and  S^1  may  be  given  using 
the  projection  tensor  A*-1  defined  by 


2.3.12 


A1 


pv 


=  g^v  +  UV 


It  is  a  simple  matter  to  check  that 


has  the  properties 
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Aaf3  U  =  0 


a 


2.3.13 


’  Aa%  - 


=  A 


a(3 


A  AaP  =  3 

ap 


AaP  A  A  >  0 
a  p  — 


where  A^  is  any  vector.  The  invariant  expressions  for  \i  ,  ,  S^V 

are  then  given  by 


2.3.14 


(i  =  T 


,HV 


=  T 


,PV 

U 

U 

p 

V 

,0,3 

A^ 

U0 

a 

3 

,0,3 

A^ 

A 

a 


The  correctness  of  equations  2.3.14  is  easily  checked  by  observing 
that  in  the  coordinate  system  where  Ua  =  (0,0,0,l)  they  reduce  to  the 
expressions  2.3.11.  Since  the  quantities  involved  are  tensors  the  equations 
are  correct  in  all  coordinate  systems. 


Next  let  us  express  M^1  in  terms  of  Ua  .  Define  p  and  by 

2.3.15  ■  p/  t  i"  j  ^  u  =  0  . 

r4. 

It  is  clear  that 


2.3.16  p  -  -  M11  U 

H* 

M  * 

and  in  the  coordinate  system  where  IT  =  (000  1)  this  becomes 

*  vr4 

2.3.17  p  =  M 

Thus,  for  the  chosen  observer,  p  is  the  number  of  particles  per  c.c. 
multiplied  by  the  rest  mass  of  one  particle.  In  the  special  coordinate 


« 
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system  =  (0, 0,0,1) 


2.3.18 


*i  *  i  „ii  * 

i  =  M  ,.  4=0 


Thus  4  represents  the  particle  flux. 


We  have  now  defined  p  and  p  ,  but  what  about  P  ?  In  equilibrium, 
spatial  isotropy  and  absence  of  viscosity  enable  us  to  define  the  thermodynamic 

pressure  as  equal  to  any  of  the  three  identical  measurable  stress  components 

12  3 

.  Away  from  equilibrium,  the  concept  of  thermodynamic  pressure 

becomes  nebulous,  since  we  are  attempting  to  apply  an  idea  of  equilibrium 

12  3 

thermodynamics  to  a  non-equilibrium  situation.  No  longer  are 
all  equal,  since  in  a  non-equilibrium  situation  bulk  viscosity  terms  appear 
in  these  components.  It  seems  that  it  is  impossible  to  separate  reversible 
(Pressure)  effects  from  irreversible  (viscosity)  effects  except  at  or  near 
equilibrium.  We  will  therefore  adopt  a  new  appraoch:  At  equilibrium  we 
are  supplied  with  a  functional  form  (The  equation  of  state)  for  P^  : 


2.3.19 


00 


=  f(u  ,  p  ) 
\  “qq  ’  f-'co  ' 


Let  us  use  this  form  to  define  P  in  the  non-equilibrium  cases 


2.3.20  P  =  f(p,  p) 

where  f  is  the  same  functional  form  as  in  2.3.19*  Since  2.2.5  only  depends 
on  the  functional  form  of  P  in  terms  of  \i  ,  p  we  will  retain  the  Gibbs 
relation 

2.3.21  T(|i,p)  dS (p,p)  =  d(^)  +  Pd(~)  ;  T  >  0  . 

Note  that  p,  p,  P,  T,  S,  depend  on  Ua  .  This  means  that  obervers 
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in  different  networks  etc.  will  observe  different  values 

|i^  ,p^  ,P^  ,T^  ,S^  and  p.^,  ?2>  ^2’  $ 2  at  t^ie  same  event  i-n  space-time. 

Only  as  t  — >  00  will  the  sets  of  observations  tend  to  agreement.  It  is 

worth  noting  here  that  the  quantities  S  ,  ^  appearing  in  2.3.21  refer 

P 

to  entropy  per  gram  of  rest  mass  and  energy  per  gram  of  rest  mass  £t  a  point. 
If  we  wish  to  express  the  relation  in  terms  of  the  total  entropy  and 

total  energy  of  a  (small)  finite  region  of  volume  ,  containing  rest 

mass  ,  we  would  proceed  as  follows:  By  2.3.21 


2.3.22 


s  i  U.  V 

Td0  =  d(^)  +  Pd(^) 


'M. 


'M, 


or 


Putting 


2.3.23 


e  = 


TdS 


+  P1V1  - 


T  M- 


=  dUx  +  PdVx  - 


thermodynamic 
potent ial/gm 


TO  dM-  . 


This  gives  the  change  in  entropy  of  a  finite  region  whose  volume  may  change, 
and  which  may  experience  gain  or  loss  of  particles  across  its  boundary.  The 
expression  ^  +  ■  -  will  appear  frequently,  so  for  future  simplification  we 
will  define  r\  ,  the  enthalpy  per  gram  by 

2.3.24  (n  +  P)  =  hP  • 

The  Gibbs  relation  then  shows 
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2.3.25 


dn  =  TdS  +  -  dP 

P 


Since 

2.3.26 


Q  = 


2  -  S 
T 


equation  2.3.25  also  implies 


2.3.27 


ae  =  i  dp  +  T,d(i) 


Having  defined  the  thermodynamic  pressure  P  we  can  now  separate 

/  UV \  uv 

the  viscous  (t  )  and  non-viscous  parts  of  the  stress  tensor  S  : 


2.3.28 


S^V  =  t^V  +  P  A^V 


Multiplication  of  2.3.28  by  and  use  of  2.3.13  and  2.3.10  give 


2.3.29 


tM-V  U  =  0  . 


UV 

We  are  now  prepared  to  express  T  in  its  final  form 


2.3.30 


uv 

T  =  T)P 


^  uv  +  pg^  +  5V  +  {V  +  T^v 


Let  us  now  investigate  the  consequences  of  the  conservation 
equations  TM'V  =  0  and  M  ,u  =  0  .  Substituting  2.3.30  into  =  0 

we  get 


2.5.31 


(nu^  puv)>v  +  PjV  g4V  +  (?V  +  £V  +  ^v)>v  -  0 


Substituting  2.3.15  into  M  ^  =  0  we  get 

f  M- 


2.3.32 


(puv)  „  -  e  „  • 
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For  brevity,  let  us  write  cp  for  Cp  UV  .  Substituting  2.3*32  into 
2.3.31  gives,  with  use  of  2.3.25 


2.3.33 


(pTS  -  r£V  +  JU11  +  £“  +  SV 

f  V  y  V  ,  V 


+  qpu^  +  P  g^v  +  pu  +  =  0 

*  *  y  s  f  v 


This  may  be  rearranged  as 


2.3.34 


TU 


([s^  +  ^  +  9  ja  .  (±)^a) 


,a  ,a- 


+  ^  +  r^V  +  ttpU^  -f-  P  A^a  +  Ca 

»  V  %  U>  y  VjU  y  vjl/ 


=  0 


Our  plan  for  the  immediate  future  is: 


1.  To  project  equation  2.3.34  on  the  time-like  vector  U 


a 


2.  To  consider  the  spacelike  projection  of  equation  2.3.34  by 


multiplying  by  A  . 

nr 


3.  To  combine  the  results  of  1  and  2  into  a  useful  expression 


for  Sa 

,a 


1.  First,  multiply  2.3*34  by  U  : 

M* 


2.3.35 


([Sff +s^a  -  (k)  f) 


T  J,a  ,a 


-  faU  -  tMV  U 

>  \±  ,v  (1 


0 
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To  simplify  this  observe  that  equations  2.3.9  and  2.3.29  imply 


2.3.36 

u  p  =  -  u4  r 

2.3.37 

u  tMV  =  -  u  t^00  , 

H  (I, a 

Also,  to  symmetrise  the  R.H.S.  of  2.3.37  observe  that  if 


2.3.38 

6  a  =  Aj  (U  +  U  ) 

cop  2  a  p  (i,v  v,p 

then 


2.3.59 

tt  M-V  ap 

U  =  -e  Q  t  . 

|i  ,  v  ap 

We  can  then  write  2.3.35  as 


2.3.40 

[stf  +  i>0  -  -  i  ^  .  etae  +  [9)>ar- 1 »/] 

2.  Now  let  us  consider  the  spacelike  portion  of  2.3.34  by  multiplying 


through  by  A^ 

Since 

2.3.41 

A^  ^  ;  AA  ;  A*  =  UA 

1 a  b  b  '  (i  *  |i  ,a  ,a 

we  get 


2.3.42 

(5“  ^  ♦  U"a  +  ^  +  P,a  ^  =  0 

To  simplify,  put 


2. 3.43 

2.3.b3 
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It  is  obvious  that 

2.3.44  6A  =  0  . 


We  now  have  the  simple  expression 

2.3.45  rpU^  =  -P  -  r)pT  5^  . 

j  cc 


To  get  this  in  shape  to  substitute  into  the  timelike  projection  2.3.40, 
multiply  it  by  l 


2.3.4 6 


u 


^°°P 

- ^  -  T6a  i 

IP  cc 


3.  We  will  now  combine  the  spacelike  2.3.46  and  timelike  projection  2.3.40 


of  T^V  =  0  : 

>V 


2.3.47 


-a  C  -  > 

T  ; 


/  cc 
(SM  + 


0, a  /(.a 


.a 


,  <!“  '  ‘  5  •«  *  *'t-  ' 


"a 


We  must  now  make  an  assumption  about  the  relation  between  the 
entropy  flux  vector  and  the  heat  flow  and  particle  flow  vectors. 

If  Sp  U^,  refer  to  the  total  entropy,  energy,  and  rest  mass 

inside  a  finite  region  of  fixed  volume  we  have 

2.3.48  TdS1  =  dU1  -  T Q  dM1 

(cf.  2.3.23).  We  will  assume  that  the  fluxes  of  entropy,  energy,  and  rest 
mass  observed  by  a  man  with  4-velocity  Ua  satisfy  the  same  linear  relation 
as  do  the  differentials  of  entropy,  energy  and  rest  mass.  i.e., 
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2. 3.49  T(S^-  SpU^)  =  -  T0(M^  -  pU^)  . 

This  may  be  rearranged  as 


2.3.50 

=  SM^  +  . 

Substituting  this  into  2.3.47  yields  our  final  equation  for  the  entropy 


production: 

2.3.51 

pH  1  ~  /„a  1  a(3  x  „ 

V  -  t  ®,a  (5  *  <  >  -  T  £apT  +  5  * 

This  is  an  important  equation  from  which  we  shall  derive  most  of  the  further 
results  of  this  chapter.  First,  observe  that  when  t  ,  T^V  ,  M*-1  , 

are  assumed  to  approach  their  perfect  fluid  forms.  Thus 


2.3.52 

’  lim  (£a  ,  f  ,TaP)  =  (0,0, C)  '  ) 

t— ♦  00 

Equation  2.3.51  then  shows  that 


2.3.53 

lim  Sa  =  0  . 

,a 

t— *  00  7 

We  were  able  to  show  (2.2.9  et  seq.)  that 


2.3.54 

8™ 
V  p 

11 

0 

so  we  now  have, 

as  expected, 

2.3.55 

lim  S  a  =  S  “ 

t-o.  ’a  00 >a 
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Section  4 1  A  Comparison  of  the  Eckai.'t  and  Landau -Li  f  shit  z  Formalisms. 

The  two  formalisms  differ  in  the  selection  of  the  network  of 
observers,  i.e.,  in  the  choice  of  Ua  .  Eckart  selects  Ua  parallel  to 
,  while  Landau  and  Lifshitz  select  Ua  parallel  to  Ea.  Thus,  in 
Ecl<  art  ’ s  rest  frame,  we  observe  no  net  particle  flow,  while  in  Landau-Lif shit z 
rest  frame,  we  observe  no  net  energy  flow  or  momentum. 

A.  The  Eckart  Formalism. 

For  Eckart  the  vector  Ua  is  defined  by 

#  =  pUa  . 


In  our  previous  work  this  implies 
S.k.2  ^  =  0  . 


The  entropy  production  is  then  given  by  2.3.51 


2.4.3 


9 

,a  t) 


1 

T  £a(3 


If  now  we  assume: 


(i)  and  Ta(3  are  linear  functions  of  (-9^-  6 p)  and  e( 

(ii)  that  the  fluid  is  isotropic 


we  are  able  to  show  (appendix  l) 


-  29  - 


2.4.4 


2.4.5 


a|3 


=  -  A. 


A7"  A^P  (U 


A.n 


U  x) 
H,  A7 


\  A 


aP 


\  >  °  a2  >  0 


.a 


=  K  £ 


aP 


(-0  -  6  > 

P 


K  >  0 


It  is  worth  noting  that  2.4.5  can  be  rewritten  by  the  use  of  2.3.45 
and  2,3.27  as 


2.4.6 


C  =  K  AaP  [(i))(J  -  Tp  ]  K  >  0  . 


This  is  the  form  of  in  Eckert' s  original  paper. 


B.  Landau-Lif shit z  Formalism. 


Landau  and  Lifshitz,  for  their  formalism  select 


2.4.7 


Ua  =  Ea 


The  first  consequence  of  this  selection  is  the  vanishing  of  ^  : 

2.4.8  T^V  Uv  =  -9  . 

Multiplication  by  U  gives,  using  2.3*7 

2.4.9  n  =  +  9 

We  then  have 

(nuV  +  pa4V  +  +  cV  +  T,iV)uv  =  -uu^ 

or  by  virtue  of  2.3*9?  2,3*29  and  2.3*7 


A<-M 


'  .  rr 
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-nir 


2.4.10 


i.e . 


r  =  o  . 


The  Entropy  production  2.3*51  then  becomes 


2.4.11 


,a 


.a 


-  V  -  T  *aP 


af3 


If  as  before  we  assume! 


(i)  and  t^p  are  linear  functions  of  Q 


a(3 


and  e 


a(3 


(ii)  the  fluid  is  isotropic 


we  are  able  to  show  (Appendix  l) 


2.4.12 


Tap  =  -A  /31  ^(U,  +  U  ,) 

1  v  A,|l  |4,V 


-  A2  AaP  U^A  ;  >  0  ;  V,  >  0 


2.4.13 


e  =  -  K  Aape 


,p 


K  >  0  . 


The  similarities  between  the  Eckart  results  and  the  Landau- 
Lifshitz  results  are  striking.  The  only  difference  appears  in  equations 
2.4.11  and  2. 4. $  ;  the  term  6^  £a  appears  in  the  former  and  not  in  the 
latter.  It  is  worthwhile  investigating  the  magnitude  of  this  term.  Let 
us  use  a  30  formalism,  and  drop  the  convention  c  =  1  .  To  be  specific,  let 
us  : 


. 


o, 
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.  use  a  coordinate  system  wherein  lja  =  (0,0,0,c) 


b.  V  =  (UX,U2,U5) 


c.  measure  p  in  gm/c.c. 


d.  measure  P  in  dynes/sq.  cm. 


e.  measure  p  in  gms/c.c 


f.  measure  T  in  K 


<,1  .  dyn-cm.  dyne 

g.  measure  r  m  — - — -  per  sc  =  — ■ L - 

sq.  cm.  cm  sec 


,  .  cm 

h.  measure  £  m  - 


2  o 

sec  K 


oTr  dyne 

i.  measure  K  m  K  — ^ — 

sec 


Immediate  results  of  these  conventions  are  I 


2.4.14 


M.  +  P/c' 


is  dimensionless; 


2.4.15 


0  -  2 


s 

7 


has  dimensions  (°K) 


A  straightforward  dimensional  analysis  of  equation  2.3 


yields 


.43  then 


o  1: 
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2.4.16 


°2  =  £div  7  +  (£' 


— »  2  — > 

V)V  +  c  div  t  + 


(V  •  V)  t 


A  similar  analysis  applied  to  2.4.8  gives 

2.4.17  f  =  K(grad  0  +  s)  . 

c 


Under  what  conditions  are  we  assured  that  grad  9  dominates 

— 2  ^  •  This  will  be  the  case  it 

c 


2.4.18 

or 


~  6  <  <  grad  9 

c 


2.4.19 


£  £  K  grad  9 


1  2 

We  know  (cf.  Taub  Phys.  Rev.  j4  ,  328  (1948))  that  F  <-  pc  for  all 

3 

’’reasonable"  fluids.  Thus 

4 

2.4.20  p  <  ip  <  r  n 

3 

rip  ~  p 

Consideration  of  2.4.24,  2.4.27,  2.4.28  shows  that  —  &  is 

c 


negligible  provided 


6 


i 
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2. Ik  22 


K  grad  0  div  V 

t  : - 

c  [i  T 


K(grad  0  *7)7 

c  (j.T 


K(V  -p)  grad  0 

~ck  [iT 


<  <  |  grad  9  | 


<  <  I  grad  G 


<  <  |  grad  0 


div  t 
c2  |iT 


<  <  |  grad  0  j 


Under  classical  conditions,  the  high  powers  of  c  in  2. 4. 22 
guarantee  satisfaction.  There  is  reason  to  believe  that  the  validity  of  the 
formalism  is  restricted  to  cases  where  2.4.22  applies.  In  these  cases 
&a  £  is  negligible  and  the  expressions  for  entropy  production  are  formally 

OL 

identical  in  Eclart  and  Landau -Lif shit z.  The  same  results  could  be  ensured 

by  restricting  ourselves  to  fluids  for  which  K,  A^,,  A^  are  sufficiently 

small.  If  we  then  retain  in  our  expression  for  S  a  only  those  terms  which 

,a 

are  linear  in  K,  Ap  A^  the  term  £a  will  be  dropped.  (it  is  clear 
from  2.3.43,  2.4.8,  2.4.7  that  5<X£a  is  quadratic  in  K,A^,  A^  .)  It  is 
obviously  a  matter  of  choice  whe^ner  we  regard  the  inequalities  2.4.22  as 
restricting  the  size  of  the  coefficients  or  as  restricting  the  permissible 
gradients . 

To  summarize  our  comparison  of  Landau-Lif shit z  and  Eclart  (as 
herein  elaborated  and  interpreted)! 


A.  Because  of  different  rest  frames  p(ECKART ) ,p(ECKART )etc . 


' 
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are  different  from  p(L-L),  p  (L-L  )  2tc  . 

B.  Both  use  the  same  functional  form  for  P  . 

C.  The  final  expressions  for  the  entropy  production  are 
formally  identical,  to  within  the  approximation  2.4.22. 

D.  The  question  of  equivalence  is  still  not  settled,  due  to  the 
(unknown  as  yet)  differences  in  p(ECKART ) ,p (L-L ),  etc.  This 
question  will  be  answered  in  the  next  section,  where  we  will 
introduce : 


Section  3-  A  New, General,  Formalism  for  Irreversible  Thermodynamics. 

We  will  leave  the  network  of  observers  arbitrary,  but  will  consider 
only  the  case  where  t  is  moderately  large  -  the  case  where  the  fluid  is 
reasonably  close  to  equilibrium.  It  is  necessary,  of  course,  to  specify 
what  we  mean  by  "moderately  large"  and  "reasonably  close".  This  presents 
our  first  difficulty. 


One  may  try  to  specify  its  closeness  by  requiring  that  the  viscous 
forces  should  be  small  compared  to  the  pressure,  e.g. 


2.5.1 


<  <  1  . 


One  also  needs  a  restriction  on  the  heat  flow;  for  example  it  may  be 
required  that  the  thermal  energy  flux  be  small  compared  to  the  mass  energy 


flux: 


► 
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— *  2  — *  3 

q  <  <  |_ic  *v  <  pc 


or 


2.5.2 


<  <  1  • 

UC5 


This  requirement  is  met  automatically  in  the  classical  limit 
where  c  -+  <»  ,  but  it  places  a  definite  restriction  on  the  range  of  validity 
of  our  relativistic  theory. 


We  however,  will  give  these  restrictions  more  geometric  form. 
Near  equilibrium,  the  vectors  Ea  and  Ma  are  near  their  equilibrium 
positions,  i.e.,  nearly  parallel.  Let 

2.5.3  Ea  -  -■  =  €a(x1  x2  X^  x^)  ;  e  =  |6a| 


€  is  dimensionless  and  we  will  take  it  as  our  standard  of  smallness.  It 
be  thought  of  as  the  pseudo-angle  between  Ea  and 


can 


We  shall  think  of  €  as  negligible  in  comparison  with  — ^ 

|1C 

1  2 

Since  (cf.  Taub)  P  <  -  pc  for  all  "reasonable"  fluids  we  have 

3 


2.5.4 


e  <  < 


pc 


<  1  . 


The  cone  obtained  by  rotating  E  and  M  about  their  bisector  will  be 

referred  to  as  the  €  cone  at  the  event  P  under  consideration. 
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Since  we  are  going  to  consider  the  fluid  only  near  equilibrium, 
we  will  require  that  the  gradient  of  ea  is  small  in  the  sense  that 
eap  / €  <  <  1  .  More  precisely,  in  analogy  with  2.4.29  we  will  assume  that 


2.5.5 


JL  ill 


b 

c (J.T 


<  <  1 


K 

T“ 

c  |iT 


.  € 


<  <  1 


Here  K  is  Fouriers  coefficient.  The  restrictions  2.5.2  insure 
that  the  €  cone  does  not  change  abruptly  in  shape  or  orientation  from 
point  to  point.  Let  us  summarise  our  restrictions: 


A.  We  allow  an  arbitrary  network  of  observers,  characterised  by  a  vector 
field  Ua(x^  x^  x^  x^)  ,  subject  to 


2.5.6 


Ua  -  ±  I  =  I  ±  e\  =  o(e) 

p  1  p 


1 


2.5.7 


i  <x  a 

U  -  E 


J  ^1  -  0(e) 


B.  The  orientation  of  Ua  relative  to  the  e  cone  should  not  vary 
abruptly  from  point  to  point.  To  wit: 


.a 


2.5.8 


K_ 

M-t  | 


<  <  1 


K_ 

HT 


1C 


a 


ic 


a 


<  <  1 


C.  The  viscosity  tensor  should  not  charge  abruptly  from  point  to  point 


2.5.9 


lT 


a3 


pTe 


<  <  1 


We  have  then,  an  infinite  set  of  possible  networks  of  observers, 
including  Eckart  and  Landau-Lifshit z,  each  characterised  by  a  vector  field 
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,a 


subject  to  2.5.3,  2.5.4,  2.5.6,  2.5.7*  Any  observer  can  decompose  T^V 
according  to  2.3.30. 


As  a  first  consequence  of  these  restrictions,  let  us  show  that  if 
the  viscous  stress  tensor  is  small  in  the  eyes  of  one  observer,  it  is  small 
for  all  observers.  To  be  specific,  assume  that  for  one  observer  (ua) 


2.5.10 


TIJV  =  OCCjP) 


a 


Consider  another  observer  with  velocity  U„  given  by 


2.5.11 


TTa  TTa 
U  =  +  A 


The  requirements  2.5-6,  2.5.7  imply 


2.5.12 


Aa  =  0(e)  . 


Now  2.5.12  implies  that 


2.5.15 


(Ua  -  l£)(U  -  U*a)  =  0(e2)  . 


Thus 


2.5.11* 


U„  u  =  -1  +  0(e  ) 

*  a 


We  then  have 


2.5.15 


C  A 

*  a 


0(€2) 


U  A  =  0(e) 
a 


2.5.16 
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The  observer  with  velocity  may  decompose  into 


a  «.a 

pU  +  5  .  We  then  have 


2.5.17 


^U*a  =  (oU»a  +  p*a  +  ?a)  U 


*a 


Since  we  have  established 


2.5.18 


P*  =  P  +  0(e  p) 


Here  we  have  used  2.5-6,  2.5  7,  2.5J2and  2.3.16. 


The  observer  with  velocity  U^a  may  also  decompose  into 

pA  +  £  •  In  view  of  2.5.15  we  tnay  write 


2.5.19 


#  -  pC  +  +  0(p«2)  , 


A  comparison  of  2.5.16  and 


2.5.20 


if 


,,a  ,a 

*  pU 


yields 


P  +  0(pc2) 


Consider  now  the  decomposition  T^V  =*  (iU^UV  +  PA^  + 
Substitution  of  2.5.11  into  this  expression  yields 


T^V  -  nU#V  +  PA£V  +^+0(6^) 


2.5.21 


+  +  (n+p^lu/  +  [C  +  (|i+P)XV]u/ 
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If  we  multiply  2.5.21  by  U*  we  get  (by  2.3.14) 

2  *5. 22  =  (i  +  0(€2|i) 

If  we  multiply  2.5.21  by  A^  we  get  (by  2.3.14) 

2-5.23  C  =  ?  +  (n  +  P)Aa  +  0(e2^)  . 

If  we  take  the  equation  of  state 
2.5.2^  P  =  P(n,p)  , 

we  may  in  view  of  2.5.22  and  2.5.18  write 
2.5*25  P*  ~  P  +  0(e  (4)  +  0(fp)  * 

It  may  be  shown  (Appendix  2)  that 

2.5.26  0(e2|i)  +  0(e2p)  *  0(e2P)  . 

Hence  we  have 


2,5.27 


p*  =  p  +  o(«2p) 


If  we  multiply  2.5.15  by  A^a  A^  we  get 


2.5.28 


a3  a3  ^  /  2_N , 

r* =  t  +  0 ( €  P) '  . 


We  have  thus  proved  that  if  the  viscosity  tensor  is  small  for  one  observer, 
it  is  small  for  any  observer.  We  have  also  proved  a  number  of  other 
signifies1-  t  results  regarding  the  differences  in  p  ,  p,  T  etc.  as 
observed  by  different  observers. 


■ 
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Let  us  collect  them  together: 


-  n+0(e^i)  ,  =  p  +  0(e2p) 


P  +  0(e2P)  ,  r£V  =  r^V  +  0(e2P) 


2.5.29 


^  +  (m-+P)^  +  0(e2n) 


e;  =  ^  +  p  ak  +  o(e"P) 


It  is  significant  that  neither  nor  ^  is  invariant  under 

changes  in  the  network  of  observers.  We  can,  however,  make  up  an  invariant 
quantity  from  the  two.  If  we  define 


2.5.30  ^  ^ 

T| 

it  is  clear  from  the  last  two  equations  in  2,5.29  that 


2.5.31 


q*  =  q^  +  o(e^p)  . 


If  we  consider  q^-  in  the  Landau-Lif shitz  frame  (U^  *  E^) 


we 


have  =  0  ,  or 


2.5.32 


* 


Thus  q^  represents  the  particle  flux  relative  to  the  energy  flux.  Since 
q^  is  independent  of  the  choice  of  U|J'  it  must  always  represent  the  particle 
flux  relative  to  the  energy  flux.  Further,  in  this  frame  Ua  =  Ea  ,  so  by 


2.5.3  and  2.3.15 
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2.5.33 


a  a 

q  *  -  pe 


It  is  worth  noting  that  2.5.32  shows  that  our  assumption  2. 3.1+9 

is  the  same  assumption  for  all  observers;  In  the  equivalent  form  2.3.50 

all  quantities  are  independent  of  U^-  .  We  may  now  write  the  expression 

for  S  a  in  its  final  form; 

,a 


2.5.34 


1 

T  ea3 


If  we  admit  the  inequalities  2.4.22  (or  if  we  think  of  K  , 
and  include  only  those  terms  which  would  be  linear  in  them), 
negligible  and  we  have 


,  as 

£>a  £  is 

a 


" small” 


2.5.35 


a 


s  =  -  e  q 


a 


,a 


a 


1 

—  €  T 

T  a,(3 


a(3 


The  requirement  that  S  a  be  non-negative,  plus  the  assumptions  of  linearity 


.a 


and  isotropy  then  yield  (See  Appendix  l) 


2.5.36 


r^v  =  -  A  Aa|J  AgV  c  „  -  A,  A(JV  AaP  €  R 

1  a(3  2  a|3 


;  \  >  o  ; 


\  >  0 


2.5.37 


qa  =  -  K  G  AaP  ;  K  >  0  . 

y  P 


The  formula  2.5.37  is  worth  commenting  upon.  In  the  classical 


case  we  have  two  laws*  Fick’s  haw  for  particle  diffusion 


i 


■ 
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2.5.38 


K 


£  = - grad  0 


and  Fouriers  law  for  heat  conduction 


2.5.39 


5  = 


K 


grad  T 


a 


In  relativity  we  have  a  single  law  q  =  -  K  0  n  A 

y  P 


a(3 


of  which 


2.5*30  and  2.5.31  are  but  different  aspects  (observed  by  different  people). 
The  relativistic  law  thus  unifies  the  two  previously  unrelated  results. 
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CHAPTER  III 

EXTENSION  OF  THE  FORMALISM  TO  A  MIXTURE  OF 
N  CHARGED  FLUIDS  IN  AN  ELECTROMAGNETIC  FIELD. 


Section  1:  Introduction. 


We  shall  assume  our  mixture  to  be  described  by  N  particle  flow 
vectors  ,  A  =  1,  2,  .  .  ,  N  ,  an  energy  momentum  tensor  T^V  and  a  skew- 

symmetric  electromagnetic  field  tensor  F^V  .  We  wj'll  assume  that  the  ratio  of 
charge  to  rest  mass  of  each  particle  of  fluid  A  is  .  If  we  select 

a  network  of  observers  by  defining  a  vector  field  Ua  we  may  write,  as  in 
Chapter  II 


3.1.1 


£ 


0  . 


Then 


3.1.2 


-  U  =  p 
A  [X  K i 


Let  us  define  the  total  particle  flux  by 


3.1.3 


*  -  >  <  ■ 


,  u  |~L  , 

If  we  define  p  and  £  by 


3.1.1* 


/  -  (>UH  +  ^  i  ^  -  0 


it  la  not  hard  to  show  that 
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3.1.5 


-M^U 


M- 


=  P 


P  = 


e  = 


The  current  flow  associated  with  fluid  A  is 


3.1.6 


=  e4  ^ 
A  A  A 


The  total  current  flow  is  then 


3.1.7 


a 


5 


A 


Let  us  define  the  electromagnetic  energy  momentum  tensor  T^^  by 


3.1.8 


^  g  f^f3v  .  1  ga3F  fpv 
(E)  pv  4  pv 


In  view  of  Maxwell's  equations (cf.  Panofsky  and  Phillips:  Classical  Electricity 
and  Magnetism,  Addison  Wesley) 


3.1.9 


T7a3  -  Ta 

f,p  •  J 


and 


3.1.10 


F  „  +F0  +  F  Q  =  0 

$7, cl  ya,  (3 


we  may  write 


3.1.11 


V 

(e),h 


-fa^j 


p 


If  we  now  define 
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3.1.12 


,nv 

(M) 


[iV  _M-V 


=  T  -  T 


(E) 


We  may  define  (i  ,  ,  s^V  by 


HUV  +  +  $V  +  SPV  =  T^} 


3.1.13 


(TU  =  0 


s^vu  =  0 

v 


As  in  Chapter  II,  invariant  expressions  for  |i  ,  ^  ,  S^V  are 


M-  =  T^v  U  U 
(M)  | jl  v 


3.1.14 


_  t°M3  Ai  TT 
C  -  -T(M)  Aa  Up 


StiV  =.  T  f  ,  AV 

(M)  a  p 


At  equilibrium  we  have  an  equation  of  state 


3.1.15 


P00  ■  f^oo'  Pi'  P2'  pn^ 

00  00  00 


As  in  Chapter  II  we  will  use  the  functional  form  f  to  define  the  thermo¬ 
dynamic  pressure  in  the  non  equilibrium  case: 


3.1.16 


PL>  P2>  •••  Pn)  * 


uv 

Having  defined  P  we  may  define  the  viscous  stress  tensor  t  by 


T^V  =  S|JV  -  P  A^V 


3.1.17 
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,pv 


We  may  now  write  in  final  form: 


3.1.18 


„PV  (  T1\1TKTV  ™  M-V  <.M-TTV  .vli  M-V 
T(M)  =  +  p)u  u  +  pg  +  Tu  +  £  Ij  +  "i  . 


Section  2:  Thermodynamic  Relations. 


At  equilibrium  we  have  the  Gibbs  relation 


3.2.1 


TdS 


=  <>  +  pd^>  -  2  T9A  d<T)  • 


Since  its  form  is  dependent  onlv  on  f  (cf.  3*1.15)  we  will  retain  it  even 
under  non-equilibrium  conditions.  If  we  define  q  and  9  by 


3.2.2 


4  = 


H  +  P 


e  =  |  -  s 

T 


we  may  show  that 


3.2.3 


d  q  =  TdS  +  -  dP  + 
P 


I  *V<T> 


and 


3.2.4 


de  ■  nd(|)  +  dP  +  V 


It  may  be  shown  that 


3.2.5 


p0  =  )  P*S 


A  A 
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Since  the  proof  is  rather  tedious  we  have  placed  it  in  appendix  3. 
Section  3 »  Conservation  equations. 

We  will  take  as  postulates 


3.3.2  S^  >  0  . 

t  M- 

Note  that  we  do  not  require  =  0  t  We  will  explicitly  include  the 

t  M’ 

possibility  of  nuclear  reactions  leading  to  a  loss  or  gain  of  rest  mass. 


In  view  of  equations  3*1 .11,  3.1.12  and  3*3*1  we  may  write 


30*5 


T 


(M),H 


J 


w 


* 


To  connect  the  entropy  flux  vector  with  the  heat  flow  and  particle  flow 
vectors  we  will  assume  as  in  Chapter  II  that  the  vectors  satisfy  the  same 
linear  relation  as  the  differentials  (3.2.1): 


T(0b8erved  entropy  flux)  =  (Observed  heat  flux) 


-  \  T0  (Observed  flux 

L->  A  of  fluid  A) 


i.e., 


3.5.4  T(S^  -  SpU^)  =  f  -  £  xeA(M^  -  pAu^) 

A 
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We  will  see  (Section  4)  that  this  assumption  does  not  depend  on 


the  choice  of  .  This  may  be  written  (cf.  3. 1.1) 


3.3.5 


-  SpU^)  = 


T0a^ 

A  A 


or 


5.3.6 


S1-1  =  SpU**  + 


eA?i  • 

A^A 


Let  us  now  consider  the  consequences  of  the  conservation  equations 


T^V  =  0  .  Using  3.3.3>  3.1.18,  3*2.1  we  may  write 


/  mo  ttV  xy  ,-V  uV  nttp  i|a  p.v  *|1  AUV 

(pTS  ,U  +  rp.  „  -  r£  v  +  £  v)U  +  £  +  t  v  +  qpU  +  P  v  A 


3.3.7 


pT  )  0A(~)  uV  + 

M  Av  p  ,  V  b  ,v 


uA 

=  J-, 


With  the  help  of  equations  3*2.5?  3*1»^?  3*1.1  and  a  bit  of  algebra  we  may 
write 


3.3.8 


y  ©A(— )  uv 

L  A  e  >v 


+  Z0A«Iv-X 

A  A 


e.  e 

A  A’v 


Substitution  of  3*3.8  into  3.3*7  yields 


-  49  - 


3.3.9 


TuMtVV  +  (f  ^>,v  +  -  C}  '  (X  9A?V\v} 


-  T(i)  jV  +  +  T)  pU^+  P  A^V  +  ^UV  +  ^ 


kT',V 


,v  ,v 


+  T  )  0  £V  u11  + 

L,  A>v  ^ 


^  - 


A 


T0  MV  UM 
A  A>V 


If  we  take  the  divergence  of  3*3.6  we  get 


3.3.10 


s 11  =  ps  uv  +  £v)  +  s(m^-^)  -  (  V  e  iv) 

,p  ,v  vTb',v  v  s  ,1-1  V/_J  As  ',v 


In  view  of  3*3.10,  equation  3.3*9  maY  be  written 


3.3.11  TU^(SVJ  -  T(i)  tV  +  S*  +  ipl?  +  P  ,  A^v  +  i 


kl 


V  VT  ,  V 


♦  tVv  ♦  T  )  eA>v  Vk  «**  +  5'tf 


V  .  I.V.J1 


F^J, 


TS,  MV  U11 
A  A'V 


We  will  next: 


(a)  Project  3*3.11  on  U 


H- 


(b)  Project  3.3*11  on  A 


|u<x 


(c)  Combine  the  results  of  (a)  and  (b)  to  get  a  simple  expression 


for  S 


a 


,a 


50 


(a)  Multiply  3.3.11  by  U  .  The  result  is 

r 


3.3.12 


-TSV„  +  T(i)  ,  tV  -  T  >  0.  e 


v  T  ,  v 


A,  v  A 


n 


+pU  +  t^V  U  =  F^JkU  + 

H  ,v  p  A  p 


TO,  M 


A  A'V 


or  in  view  of 


3.3.13 


u  ^  =  -U^ 

M.  bp 


U 

M-  ,v 


£apT 


a  3 


Cf.  (2.3.36,  2.3.37)  we  may  put 


3.3.14 


=  A)  tv 

S,U  V,v^ 


Z  •»,.« 


v  .  I  >  t  -I  T^e 

A  T  U  V  T  €A|j 


A 


i  F**\u 

T  A  p 


9  MV 
A  A’V 


(b)  Multiplying  3«3.H  by  A  we  get 

r* 


3.3.15 


/  uv  „A  :u  .A  „A1tv  ..A  t*.A  A 

(T,  v  A|a  +  4 ^  +  ?  U  „  +  U  P  >  +  7PD  +  P  v  A 


Av 


,v  ,v 


=  F^J  AA 
V  p 


Putting 


3.3.16 


„  ,A  -u.  A  ,A„v  „A  ,v  uv  aA 

nPIb  -  +  £  u,v  +u,v^  +  ,  v  Au 


we  have,  upon  multiplication  of  3-3-15  by  ^ 
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3.3.17 

x  cvp  f^vj  aa  £ 

A-V  .  v  _v„  v  p 

U  ^  -  - ■-  -  T5  Cv  +  - * -  • 

IP  IIP 

(c)  We  will  now  combine  3*3. 14  and  3«3.17  to  get  an  expression  for  S*1 


3.3.18 

1  V1  P  F^J  a\;x 

oP  A\  \  .  tv  „v  v  la^A 

V  ■  <t\v«  -  l  eA,v?A  +  "qpT  ^  +  6  -  * 

A  ^pT 

pMVj  u  1  -v 

-  -l>e„  9a  Mv  . 

T  T  v  A  A  A,v 

Section  4:  Invariance  under  changes  of  Ua  t 


As  in  Chapter  I,  let  us  define  e  by 

3-4.1 

|Ea-i#|  -  e  . 

We  will 

assume  that 

3-4.2 

3.4.3 

|ua  -  Ea|  =  |  i  ^  |  =  0(e) 

To  this  we  will  for  consistency  add  the  requirement 

3. ft.1*  FtiV  =  o(e)  ;  T4V  =  0(e)  . 
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If  this  were  not  the  case,  electrons  and  protons  would  be  accelerated 

powerfully  in  opposite  directions,  MI,  and  MI  would  differ 

’  Electrons  Protons 

by  an  appreciable  amount,  and  3*4.2  could  not  be  satisfied.  This  requirement 

2 

shows  that  to  within  0(e  ) 


3.4.5 


F^J  U 
H  v 


is  independent  of  rest  frame 


It  is  easy  to  see  that  under  a  change  of  rest  frame 


3.4.6 


TTa  TTa 

U  =  U*  4*  A  «, 


the  quantities  T,  S,  p  ,  0  ,  p.  are  invariant  within  0(e  )  while 
transform  as 


3.4.7 


^  +  (n+P)^  +  0(£2^) 
it  =  ^  +  p.  ^  +o(£2p) 


Thus,  dropping  terms  in  0(e  ) 


(SpU^1  +  ^  -  £  eAe^)#  =  SpU^  -  SpA^  +  ±£v  +  ^ 


Va 


•  I  v/“ 


or  in  view  of  3*2.2,  3*2.5 
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3.4.8 


(SpU1-1  + 


X  ®A?A>*  -  S*U"  +  ¥*  '  I  «A«A  • 


Our  assumption  about  the  form  of  the  entropy  flux  vector  is 

a 

thus  undependent  of  the  choice  of  U  .  Every  quantity  which  went  into 
the  equation  3*3.18  was  independent  of  Ua  ,  so  the  R.H.S.  of  3*3.18  must 
be  independent  of  Ua  .  It  is  illuminating  to  transform  3.3*18  into  a  form 
which  is  obviously  independent  of  Ua  . 

First,  observe  that  in  view  of  3*4.2,  the  quantity 


5.4.9 


=  fv  .  ^ 

4  A  ^A  Pr|  4 


fJ. 

is  independent  of  the  choice  of  U  .  In  the  L.L. 
have  =  0  and  3*4.3  says  that  q1^  is  the  A*"*1 
to  the  momentum. 


frame  (U^-  =  E^)  we 
particle  flow  relative 


Next  define  e 


and  by 


3.4.10  f  =  peU^  +  /  ,  V'^U  =  0  . 

H' 

By  considering  a  frame  wherein  =  (Q,G,Q,l)  we  see  that  ^  Is 
observed  current.  It  is  easy  to  see  that 


3.4.11 


ep  * 


eApA 
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Under  a  transformation  U^  =  U^  4-  A^  we  have 


C  =  /  +  peA' 


3.4.12 


p  M’ 


=  ^  +  (n+p)^ 


Thus  the  quantity 


3.4.13 


,|a  ,u  e  „li 

1  =  r  -  -  r 


eA(qA} 


is  independent  of  (to  within  0€^)  .  In  the  L-L  frame  (U^  =  E^)  we 

have  =  0  ,  so  is  the  current  relative  to  the  momentum. 


The  electro  magnetic  terms  of  3«3»l8  are 


1  XVj  ^  h 

if  v  u  A 

”t  V  t  p 


+  f“vJ  U 
v  u 


In  view  of  3.4.4,  3.4.5,  3*4.7  and  the  skew  symmetry  of  F^V  we 


have 


3.4.14 


,  /  ^vj  aA  U 

JL  (  V  [l  A 

TV  TP 


+  F^V  J  U 

v  H, 


|  ( fm'vj  (— t;  +  u  ) 

T  V  Vvqpbp  [X 


=(FtiVjvU  )  +  no 

1  V  (J.  Tjp 


k^k  A 


V  • 


The  second  term  is  quadratic  in  the  small  fluxes  *  We  wil1  therefore 


put 
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3.4,15 


AP  =  6U 


1  xv  Eek 


fF 


TIP 


and  rewrite  3.3*18  as 


3.4.16 


li 


‘  ^  0A,V^A 


2  (-‘'l 


A  r|pT 


.  A^t  -  Ir^e 
V  A(i 


0  MV 
A  A'V 


I  FpVj  u 

T  v  |i 


We  may  rewrite  the  quantity  in  square  brackets  in  terms  of  q  : 

x\ 


(-)  tV  -  \  w  . 

VT;,v^>  ^  A,  v  A  qpT 

A 


..I!  +  7^  ^ 


,PA.,  \ 


?  ((t\v  +  +  i,  X  9Al'p'>,J 


3.4.17 


V?I 


+iv- 

1 


0 


A,v  P 


Vv  e; 


SA  v  (?A  -  ~  cv) 

'  1 


0.  q, 
A,v  A 
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Here  we  have  used  equations  3.2.5,  3.2.4,  3.4.3.  We  may  now  write  S 


as 


3.4.18 


M- 


0.  q^ 

A,  v  A 


1  An 
t  T  €An 


0  MV 
A  A,v 


4  A  J  U,  +  A1^ 

t  n  a 


If  we  assume  the  coefficients  of  viscosity,  heat  conduction,  electrical 
conductivity  and  diffusion  as  small,  we  may  drop  A^£  .  In  view  of  3*4.7 

we  then  have 


3.4.19 


■  M- 


=  -  |  \  T0.  MV  -  |  \  (T0  +  e  )  q 

T  L  A  A>v  T  L  A’v  A  v  A 


v 

A 


1  A(j. 

-  —  TP 

t  An  ‘ 


The  terms  in  3*4.13  may  be  interpreted  as  follows: 


(a) 


0  M  is  the  entropy  production  due  to  the  chemical 
A  A’V 


reaction. 


_A*__  x  v 


(b)  (T9S  „  +  e^F”  U-^)q^  is  the  entropy  production  due  to 


L  A>v 

A 


diffusion  (diffusion  of  charged  components  =  current  flow) 


/  \  1  An 

-  T  T  % 


is  the  entropy  production  due  to  viscous  forces. 
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If  we  assume  a  single  reaction  we  may  write 


3.4.20 


M  =  v .  V  , 

1  j  A  " 


where  V  is  the  velocity  of  the  reaction.  Putting 


3.4.21 


9  =  )  T9ava 


we  have 


3.4.22 


T0a  M 

A  A,v 


=  <pv 


To  separate  bulk  viscosity  effects  from  ordinary  viscosity  we  may  put 


3.4.23 


1  a 
3  Ta 


;  «  Tap  -  ^ 


,a 


The  expression  for  S  may  then  be  written 


a 


3.4.24 


1  =  — 
,[i  ~  T 


**  +  I  (T\v  +  +  +  «eaj 


The  terms  in  are  of  3  kinds:  products  of  scalars,  products 

1 

of  vectors  and  products  of  tensors.  If  we  assume  that  the  system  is 
isotropic,  the  terms  of  different  species  must  be  independent. 


Section  5:  Conclusions  and  Interpretation. 


A.  Let  us  assume  that  V  and  jt  are  linearly  dependent  on  9  and  € 


a 

a 


Then 
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V  =  acp  +  be 


a 

a 


5.5.1 


a 

rt  =  c  €  +  dcp 

a 


The  contribution  of  these  scalar  quantities  to  the  entropy 


production  is  then 


5.5.2 


1  ,  2  ,  a  ,  a  ,  cu2  , 

-  -  {acp  +  be^cp  +  de^cp  +  c(ej  }  . 


a  a 


a 


The  requirement  that  >0  then  implies  that 


2  2 

-  (aa  +  (b+d)a3  +  c3  } 


be  positive  definite. 


B.  Let  us  assume  that  q  is  a  linear  function  of  T 0n  +  e^F^'th 

3,v  3  v  A 

Then  we  have 


5.5.3 


\  (L  )  (T@  +  e_FA*UO 

/  AB  B ,  v  B  V  A 

B 


Let  us  introduce  the  diffusion  flux  D^T  of  the  fluid  A  relative 


to  the  total  mass  flux  : 


3.5.4 


Da  =  < 

A  A 


A 


Since,  by  definition 


’A 


_A  p. 

)  ^ 

n 


3.5.5 
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we  may  put 


5.5.6 


D.  =  q 
A  nA 


r  • 


Summation  of  5*5. 5  yields 


3.5.7 


■  e11 -7  5“ 


We  may  then  combine  3.5.3,  3*5.6  and  3.5.7  to  get 


3.5*8 


D 


Jl 

AS 


B 


Vnf  +epFX*UTv 

CBJ  B,p  B  p.  A 


If  we  define  NA„  by 

AB 


3.5.9 


pA 

N  =  L  -  — 

AB  AB  p 


CB 


we  get 

3.5.10 


D 


M 


=  \  nat>(ts  +  epFA#u\) 
Lx  45  B,P  B  [i  V 


If  no  magnetic  field  is  present  we  get 


3.5.H 


)/aBT9B,(J 


which  is  the  relativistic  generalization  of  Fick's  law  of  diffusion,  including 
the  possibility  of  cross  effect  between  fluids  of  different  kinds. 
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Next  let  us  see  what  3*5*3  implies  about  the  current  vector. 
If  we  take  the  vector  U  =  Ea  ,  j^1  is  the  observed  current.  We 

have 


3*5*12 


P 


(i 

e  .q. 
AnA 


=  )  e  )  La_(T0_  +  e/'lK ) 

A  AB  B,(i  B  p  A' 


A  B 


5*5*13 


P 


+  CTFm*U^ 
A  AB  B,|_t  P  A 


A  B 


where 


3.5.14 


=  L  e*L  lab  eB  • 

A  B 


For  a  mixture  in  chemical  equilibrium 
the  relativistic  Ohms  law 


0)  we  have 


3*5*15 


A. 

j  =  o-F  IK 
P  p  A 


In  a  mixture  not  in  chemical  equilibrium  3*5*13  shows  that  the 
current  depends  on  the  gradients  of  the  chemical  potentials  as  well  as  the 
electric  field. 


C.  Let  us  assume  that 


-a|3 

T 


is  a  linear  function  of  €  _ 

<x[3 


Isotropy  then 


leads  us  to 
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3.5.16 


iap  =  nA^A^fU  .  +  W.  )  -|nAap  U7 

1  7,5  b,7J  3  ,7 


In  view  of  3.4.17,  3.5*1  we  may  put 


3*5.17 


a3 


=  n  Aa7  Af3b{U  +  U  }  +  (C  -  fn)  Aaf3  U7  +  depU7 


7,6  6,7 


,7 


The  first  two  terms  are  the  usual  terms  in  the  stress-tensor;  the 
last  one  expresses  the  possibility  of  cross  effects  between  the  bulk  viscosity 
and  the  chemical  reaction.  We  may  note  here  that  similar  work  (with  special 
selection  of  the  4  velocity  Ua)  on  the  relativistic  irreversible  thermo¬ 
dynamics  of  mixtures  has  been  done  by  Stueckelberg  [2]  and  Kluitenberg, 

De  Groot  and  Mazur  [6],  [9]>  [10]. 
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CHAPTER  IV 

THE  STRUCTURE  OF  PLANE  SHOCK  WAVES 

Section  1:  Introduction. 

A  shock  wave  is  a  physical  phenomenon  which  may  be  produced  by  an 
apparatus  as  shown  in  the  diagram: 


When  the  valve  is  opened,  fluid  in  the  tank  rushes  into  the  pipe, 
causing  a  phenomena  called  a  shock  wave: 
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If  the  pipe  is  long  enough  there  is  a  uniform  region  at  (high) 
pressure  ,  entropy  ,  velocity  ;  a  transition  layer  or  shock 

wave  in  which  P,  S,  U  vary;  and  a  uniform  region  (which  the  shock  wave 
has  not  yet  reached)  at  pressure  P^  ,  entropy  ,  velocity  . 

A  relativistic  treatment  of  shock  waves  is  of  current  interest. 
In  the  theory  of  gravitational  collapse  (recently  proposed  by  Hoyle  and 
Fowler  to  account  for  quasi-stellar  radio  sources)  the  shocks  generated 
by  novae  and  supernovae  should  certainly  be  described  relativistically . 


Our  purpose  is  to  investigate  the  variation  of  P  and  U  across 
the  shock  wave,.  To  make  the  mathematics  practicable  we  will  consider  an 
infinite  plane  shock  wave  in  an  otherwise  uniform  medium.  We  adopt  a 
coordinate  system  wherein  the  shock  layer  is  at  rest.  Assuming  the  plane 
of  the  shock  wave  perpendicular  to  the  X  axis,  we  will  have  all  quantities 
as  function  of  X  alone. 


y^e^ 

^  h  o  <-  ^ 

V/  (A  V  c 

. JA _ 

P  €  s  5  u  r  e 

Pm 

tr\r<o  fi  uj 
S  (X) 

y  e  I  o  c.  I  !  * 


r  <2  cj  i  a  n 

&  +  C  O  n  $  q 

I  o  u;  i>  rc  Jsc 

i~  j'op  cj 
, 

2  C?yi  s 
V  £.  I  o  c  )  4  cj 
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The  observer  runs  along  with  the  shock  wave,  so  that  the  quantities 


will  not  depend  on  time. 


Section  2;  The  Conservation  Equations. 


For  one  dimensional  flow  (steady  in  the  frame  wherein  the  shock 
wave  is  at  rest)  we  have  (x”*  =  x) 


4.2.1 


!M-=  o 

cx 


6t 


pi 

-  =  0 


)x 


with  first  integrals 


4.2.2 


m1  =  J  ;  T^1  = 


where  J  and  are  constant.  Writing  these  out  explicitly 


4.2.3 


pU1  +  =  J 


4.2.4 


rjplAl1  +  Pg^1  +  U^V  +  U1^  +  t^1  = 


The  obvious  choice  of  a  is  that  which  makes  £  =  0  .  Thus,  select 

UM  =  i  .  For  this  choice  of  ,  as  we  have  shown  in  Chapter  II 

P 


4.2.5 


p  =-l  ^P  =  _K0  AiP  .  ^  =  o 


4.2.6 


T^V  =  -  A  A^  A^V  e  -  A  A  A  •  € 

1  a3  2  a3 


-  X  A^V  AaP, 


For  one  dimensional  flow  we  have 
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u1  =  U  U3  =  0 

4.2.7  — 

U2  =  0  u4  =  7(1  +  u2) 

Substituting  these  into  4.2.5  and  4.2.6  then  yields 

b.2.8  c1  =  nK(i+u2)||  ;  c,k  =  nKu/(T+ij^)|| 

t11  -  -(VV(uu2)  I 

4.2.9 

TU1  =  -(A  +A  )u  Jut j2 

12  dx 

If  we  simplify  by  putting 


4.2.10  a(x)  =  T|K  b(x)  = 

we  get,  upon  substitution  into  4.2.3  and  4.2.4  the  three  equations: 

4.2.11  pU  =  J 

4.2.12  P  +  qpU2  +  2aU(l+U2)  ~  -  b(l+U2)  =  K1 

4.2.13  r)pU  7T+U2  +  a(l+2U2)  \l  1+U2  -  bU  7l+U2  =  K2 


These  equations  must  hold  in  the  undisturbed  region  1  ahead  of 

6u 

the  shockwave,  wherein  =  —  =  0  .  Hence,  denoting  quantities 

OX  cx 
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measured  in  this  region  by  1  we  have 


4.2.14 


J  -  biui 


4.2.15 


K1=  HjPjU2 


4.2.16 


k2= 


v,1ur^*i 


The  final  form  of  the  equations  to  be  solved  for  P(x)  and  S(x)  is 
then 


4.2.17 


pU  =  J 


4.2.18 


-p  +  -  -1)j2  +  2au(l+u^)  -  b(l+U*) 

1  p  p^  dx  ' 


2n  80  ,,  „2n  8u 

8x 


=  0 


4.2.19 


l+u2  -  TJ1  n/1+U2  +  a(l+2U2)  Vl+u2  ||  -  bu  Vl+U2  ||  =  0 


Section  3«  The  Relativistic  Rankinc-Hugoniot  Equations . 


The  equations  4.2.17  -  19  hold  in  the  portion  of  the  fluid  behind 

the  shockwave  (2)  .  In  this  region  =  0  .  Thus,  denoting 

v  '  dx  dx  70 

the  quantities  measured  in  region  2  by  a  subscript  2  we  have 


4.3.20 


P1U1  =  P2U2 


4.3.21 


(—  -  ^i)j2  +  (P2  -  Pt)  =  0 


4.3.22 


4  ..2  /.  ..2 


+  Up  -  vl  +  U^  -  0 
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no  \ 

Multiply  4.3*21  by  —  +  — 

p2  Px 

and  subtract.  The  result  is 


multiply  4.3.22  by 


+ 


M.23  ({  -  {)  -  (p2  -  p^  +  ^j  =  0  . 

The  equatior s  4- 3*20,  4.3*21  and  4.3.23  are  the  Relativistic  Rankin-Hugoniot 

equations  (cf.  Taub  [4]).  The  equation  4.3.23  implies  an  interesting 

relation  between  the  pressure  differential  -  P^  and  the  entropy 

2  2  n  ^1 

differential  -  S,  .  We  will  expand  n  -  n,  and  -1  +  —  in  powers 

2  1  1  '1  p 

of  S  -  and  P  -  P^  about  a  point  in  the  region  1.  Noting  the 
thermodynamic  relations 


4.3.24 


( 


$  ri/p 

dP 


( 


we  may  write 


u.5-25 


Tl-T)!  =  T^S-S^  +  V^P-Pj) 


(P-P^2 

21 


( 


(P-P,)5  k 

y  ~  +0  (S-S^  +0  (P-P^ 


+ 
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4.3.2 6 


H 

P 


2  —  + 
P1 


LV*'.(Sypj 

1 


(s-s^ 


[vf +  ( I  )s  ](r  -pl) 


4^) 

\  oP  / 


(p-p^' 


dp  /s]  2 : 


H  ofS-Sj)  +  o(P-P1)3  . 


Note  that  in  4.3*2  5  and  4.3*26,  o  (S-  S,.)  is  a  sum  of  terms  of  form 
A(S-S1)m(P-P1)n  ;  m,n  >  1  . 


We  then  have 


4.3.27 


-  (p-px)(  *  +  ^  )  -  z^is-sj  -  (  ^ 


(p-p^ 


3! 


4 


+  °(S-S1 )  +  ofP-P^*  . 


In  view  of  equation  4.3.23,  we  must  conclude  that 


4.3.28 


S2'S1  12 


(p  p)3 

Vl  \  /s,  ^  x 


i.e.,  the  total  entropy  discontinuity  is  of  the  third  order  in  the  total 
pressure  discontinuity. 

Section  4:  Approximation  to  the  differential  equations. 


In  this  section  we  will  derive  an  approximation  for  the  equations 
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4.2.17  -  4.2.19  in  a  form  suitable  for  solution.  Let  us  define  e  by 


4.4.1 


Let  us  assume  that  the  major  part  of  the  variation  in  P  takes  place  over 
a  distance  of  order  L  .  We  shall  assume  that 


4.4.2  e  -*■  0  K-  >  L  “+°o  . 

1  2  1  e 

For  out  approximation  we  shall  retain  terms  of  order  e  ,  “  ,  e  ,  ,  “ 

L 

and  drop  all  higher  ones.  We  then  have,  assuming  that  P  varies 
smoothly  and  monotonically  through  the  shock  wave 


P  -  P, 


0(e) 


4.4.3 


f.  ^  =  O(^)  . 


Assuming  the  variation  of  S  inside 


and  subject  to  only  a  few  maxima  and  minima  we 


the  shockwave  is  smooth 
have 


4.4.4 


dS 

dx 


max 


S-S, 


1  'max 


Let  us  now  begin  by  multiplying  4.2.10  by  r/p 

/ - p  / — “p 

by  (t)s/1+U  +  r)1  Vl+Up/J  ;  and  subtracting. 


+  rij/p^  i  multiplying  4.2,19 
The  result  is 


4.4.5  T,2-^  -  (P-P^fVP  +  ypx)  +  2(ni/j)(l+uf)ai  §  +  o(S-S^)  =  0  . 
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In  view  of  4.3.27  we  have 


4.4.6 


-(l+Uj) 

S-S1  +  o(S-S^ )  =  ~ 


be  \  dp 
5P  Js  dx 


/  60  \  dS 
\  6s  J  dx 

1 


Using  4.4.4  we  have 


4.4.7 


s-s 


1  'max 


(const  j  J 


dP  j 

dxl nax 


0( e/L) 


Thus,  by  4.4.4  we  have 


4.4.8 


and 


dS 

dx 


To  our  degree  of  approximation,  we  may  therefore  drop  o(S-S^) 
Equation  4.4.6  then  becomes 


4.4.9 


Next  let  us  approximate  4.2.18.  Using  4.3.26  we  get,  to  the 


desired  degree  of  approximation, 
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4.4.10  f 


6s 


,f-v  *  (  ^  )s<-,> .  ( 


.2  ,  x  (P-P1)i 

2 _ nZeL.  \ : _ -, 


6F 


A,  2’ 


+  (P-P^  +  2*^  ^  ^  ^  ^  (I+U2) 


2X  d? 
dx 


=  0 


Note  that  we  have  used  the  relation 


4.4.11 


dU  _  dV  _  /  6\  \  dP 

d£  =  J  aJ  =  J  (  S?  l  dl 


To  get  the  equations  in  a  soluble  form  let  us  now  substitute  4.4.5  into 
4  4,6.  The  result  is 


4,4.12 


dP 

dx 


6s  J?l  Tl  \  //s1 


1  !\  Jsl  Jsl 


-  J 


1+U 


,2 


^  j 


(F-P,)- 


sLiZal  22£ 


6p‘ 


2  A  2  J 


Thir  may  be  simplified  by  observing  that  the  right  hand  side  is  a  polynomial 
of  degree  2  in  P  with  the  sseros  and  (left  hand  side  vanishes 

in  region  1  and  region  2).  Tims 


4.4.13 


dP 

dx 


-J2 

2 


dln4 


dP 


2  is. 


(p-p,)(p-pJ 


(l+U?) 


Vp)n-  JL  Z'  x', 

.‘4*®  jP  t^Sp^ 


'  b9  \ 
■iui  i  3p  ;c 


J  +  2a, U,  f  -r-  )  -  Jb, 


/d/ 

i(  6? 


/ 


s. 
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In  view  of  the  thermodynamic  identities 


(  d(n/p)  A 

\  dS  \ 


V  T 
11 


1 


/  be 
\  5P 


dv 

c)S 


1 


this  may  be  rewritten  as 


it. 4. 14  ^ 

dx 


J 

2 


-^^\(p-Pi)(p-p2) 
ap  /si 


(1  +  up 


alTl 


\ 


be 

Ws, 


h  f  3v 

bi(  5p 


It  is  interesting  to  compare  this  with  the  classical  expression  for  -r* 

d  X 

as  given  by  Landau-Lif shit z  [3] 


dP 

dx 


J 

2 


(p-p1)(p-p2) 


The  expressions  are  very  similar 

of  V  by  and  of  T  by  0 

P 

already  been  commented  on  (cf. 
classical  approximation  equal  to 


,  the  only  difference  being  the  replacement 


.  The  replacement  of  V  by  (-2)  has 

P 

be 

Israel  [5])*  The  term  (^p)g  is  the 


-  —  (^) 

„2  vdP's, 

T  1 

il 


as  0  =  ^  -  -p  ,  in 
C 


tv  • 
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the  classical  limit  n  1  and  c  co  we  have 


S  =  (£) 


Thus 


a  TfM 

1  ll  dp 


alTl(  2 

\  T 


dT 

dP 


ai  /  aT\2 


=  I< 


/  dT 

\  Sp 


Thus  our  equation  4,4.10  reduces,  as  it  should,  to  the  classical  expression 
in  the  classical  limit. 

Section  5t  Conclusion, 

To  continue,  in  equation  4,4.10  put 


4.5.1 


j  /  a...  la/pj, 

2V  sp2 


Js 


^  +  [v  ( I 


\2 


•  <  8 ) 


S,  j 


We  then  have 


4,5.2 


dP 

dx 


A(P  -  P1)(P  -  P2) 
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or 

4.5.3 


_ dP 

(P-P1)(P-P2) 


Adx 


Integrating,  and  putting 


Ax  = 


=  0  where  P  =  ■|-(P2  +  P^) 


1  ,  -1 

lTp  '  p'  ^  tanh  - 

2(P2  V  MW 


we  have 


or, 


b.5.b  P  =  ^(P£  +  Pt)  +  i(P2  -  Pt)  tanh^Pg-P^Ax) 

The  consequences  of  4.5*4  are  best  seen  in  a  graph. 


V  -> 


)rlo»J  (j'i  -  1 


-  75  - 


As  the  diagram  shows,  the  higher  the  pressure  differential,  the  more 
sharply  defined  and  narrow  is  the  region  which  could  be  called  a  "shock 
wave" . 


The  equations  4.5.*+  and  4.5.6  may  be.  combined  to  get  an  expression 
for  S  -  : 


4.5.5 


sech2  (|(P2-P1)Ax) 


The  graph  of  S  is 


A 

s 


H  i  cj  b  P  r  c  i  s  c/  r*  & 
D  i  U  c*  r  £  n  c  & 

Low  PrC 


X 
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As  the  graph  shows  the  entropy  has  a  maximum  inside  the  shock  wave. 
4.5.5  implies  that  =  0  .  The  reason  for  this,  of  course,  is  that 

is  of  third  order  in  the  pressure  differential,  while  we  have 
only  retained  terms  of  order  2  in  the  pressure  differential.  Hence,  to 
our  degree  of  approximation,  is  in  fact  0. 

The  length  L  over  which  the  majority  of  the  change  in  pressure 
takes  place  is,  according  to  4.5.1  and  4.5.5>  proportional  to  [i(P2~P^)A]  * 
The  A  is  a  constant  depending  only  on  the  thermodynamic  properties  of  the 
undisturbed  fluid.  Hence  the  length  L  is  proportional  to  -•  ,  and  our 

approximation  procedure  is  seen  to  be  merely  that  of  ignoring  quantities  of 

3 

order  e 


r; 
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CHAPTER  V 

TOLMAN ’ S  THEOREM 

Section  1:  Introduction. 

In  this  chapter  we  will  consider  the  following  problem:  in  a 
particular  coordinate  system  (*)  the  gravitational  field  appears  static. 
What  may  be  said  about  the  temperature?  Tolmai  [6]  has  shown  that  if  thermal 
equilibrium  is  assumed  then 

5.1.1  *s/gi  ^  T  =  constant  . 

We  will  show  that  this  assumption  is  unnecessary.  All  that  is  needed  is  the 
reasonable  assumption  that  as  t  ->  «>  all  heat  flow  ceases.  We  are  then 
enabled  to  prove  that  thermal  equilibrium  is  implied  by  the  fact  that  the 
field  is  static. 

Section  2 i  The  Form  of  the  Energy  Momentum  Tensor. 

In  the  starred  coordinate  system,  the  field  is  static.  Thus 

5.2.1  '  8  0  )  *  0  , 

dx 

Define  a  vector  Ua  by 

5.2.2  U®  *  (0,0,0,  — i—  )  . 

'I'ikb 


It  is  obvious  that 
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5.2.3 


iAj  =  -l 

a 


We  then  have  the  result 


5.2.4 


.a 


u  U  =  0  . 

>P  a 


Also 


5.2.5 


♦a  a  B 

U  U  =  U  0UPU  =  0  . 

a  ,3  a 


Observe  that 


5.2.6 


dU 


u  Q  +  u  UQ 

a,  3  a  3 


6x 


a 

3 


au 


L°  J  u  +  ! — ~ 

ap  0  |^XP 


-  }  u  1  upu. 

Q0  |ij  3 


or 


5.2.7 


*  8u  ,  r8u 

u  +  u  u  =  — ^  ^  )  u,  + 

a, 3  a  3  ^3  la3J  ’■ 


t  -  O  \ 


4 

u  u. 


In  view  of  5*2.1  this  gives 


5.2.7 


U  Q  +  U  U.  =  0 

a,  3  a  3 


Contraction  of  5.2.7  gives,  with  5.2.5 


5.2.8 


ua  =  0 

,a 


Consider  the  fact  that 


U  rMa 
(i*a37 


Ua,37  ”  Uop,73 


5.2.9 
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Multiplication  by  ga3  gives 


5.2.10 


U 

M-  7 


U3  -  U3 
,37  ,73 


*  u 


,73 


by  5.2.8 


=  - (U3U  ) 


7',  3 


by  5.2.7 


-  -(6VUr 


by  5.2.6  and  5.2.5  . 


We  may  then  write 


5.2.11 


U  [r|J  -  ^  R] 


H  7 


-  iR)°. 


In  view  of  Einstein’s  equation 


5.2.18  xT^  +  [R®3  -  R]  -  0 

5.2,11  enables  us  to  write 


5.2.13  fy.  .  AUa  . 

Thus  w®  have  the  result  that  the  velocity  Ua  of  an  observer  who  seas  the 
gravitational  field  as  static  is  the  timelike  eigenvector  of  the  energy 
momentum  tensor,  According  to  equations  2,3.8  -  10  and  2,3.14,  2.3.28  we 


have 
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5.2.14 


TfiV  =  nuV  +  P  A^V  +  jV  +  £  V  +  T^V 


where 


\i  =  T^VU  U 
M-  v 


C“  -  -TaP  A^  Ue 
a  p 


T^V  =  Ta|3  ^  A*  -  P  A“V 
a  p 


Since  U*1  is  the  velocity  of  the  observer  in  the  Landau-Lif shit z 

frame  (U^1  =  =  timelike  eigenvector  of  T^V)  we  may  use  equation  2.4.10 

and  conclude  that 


5.2.15 


r  =  0  . 


If  we  use  equation  2.4.4  in  conjunction  with  5.2.4,  5*2.5  we  get 


5.2.16 


t^V  =  0  . 


The  energy  momentum  tensor  is  now 


5.2.17 


T^V  =  (|i+P)UKlT  +  Pg 


MV 


We  know  that 


8T 


,a3  * 


0 


dx 


in  view  of  5.2.12  and  5.2.1.  Then  5.2.17  and  the  conservation  equation 

TaP 


0  give 
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5.2.18  ^  *  5P  o  . 

dx  dx  dx 

Since  p  =  p(P,|i)  we  have 

5.2.19  ^  *  0  . 

dx 

Expressed  covariantly  5*2.18  and  5.2.19  are 


5*2.20 


P  U*1  =  [i 


The  conservation  equation 


0  then  gives 


5.2.21 


(|i+P)U  +  P  =  0  . 

H  ,|i 


Since 


0  *  U  Ua 

[x  ji,a 


U 


a 


# 

w 


1 

84U 


we  may  write 


5.2.22 
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Equation  5.2.21  then  gives  us 


5.2.23 


dP 

^x1 


1  *  S  iWIgT j 


l-l  |i+P 


8xK 


From  the  thermodynamic  relation  2.3.27,  we  have 


5.2.24 


1  d0  =  J-  dP  -  d(in  T)  . 

T)  |i+P  '  ' 


Combining  the  two  above  equations  gives 


5.2.25 


-  0  "  "  '  in^"g).).  T) 


Section  3»  Proof  of  Tolir.an's  Theorem. 

Multiplication  of  5.2.24  by  U  gives,  with  5*2.20,  the  result 

f* 

5.3.1  ( t)  =  0  . 

Now  we  know  that  V-g.^  T  Is  constant  along  the  world  lines  of 
the  observer  who  see  the  field  as  static.  If  we  can  show  that  at  some 
time  \T -g^  T  is  the  same  for  all  observers,  we  know  that  it  is  constant 
throughout  space  time.  To  do  this  we  will  make  the  assumption  (reasonable 
in  view  of  the  second  law  of  thermodynamics)  that  at  t  =  00  thermal 
equilibrium  is  reached.  Thus  at  t  =  00  f 

e 


5.3.2 


0 


,  -  •  -  7,  *  .  5s 


-  83  - 


Equation  5.2.24  and  5.3*1  then  show  that  sA-g^  T  is  universally  constant. 

Further  5.2.24  assures  us  that  0  is  universally  constant.  Since  the  heat 

flow  is  dependent  on  9  ,  we  have,  universally,  no  heat  flow.  The  system, 

y 

then,  must  be  in  thermal  equilibrium.  Since 

Q 

5.3.3  9  =  —  =  constant  (G  =  Gibb’s  Function) 

we  have 

5.3.4  G  s/-g1^  =  constant 
which  is  Klein’s  theorem,  [7] 

To  sum  up:  we  have  shown  that  the  assumption  of  a  static  field, 
plus  the  assumption  of  thermal  equilibrium  at  t  ■  »  are  sufficient  to 
show 

1,  nT-844  T  -  constant  (Tolman’s  Theorem) 

2.  A  44  G  *  constant  (Klein's  Theorem) 


3* 


Thermal  equilibrium  (No  Heat  Flow)  at  all  timee. 


. 

,  '  •  >  • 
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Al.l 

Let  us 
Then 

A1.2 

A1.3 

Since 

A1 .14- 

Take  a 

A1.5 

A1 .6 


APPENDIX  1 


THE  FORM  OF  THE  "VISCOUS  STRESS  TENSOR 


We  have  the  result 


A  .  a  1  af3  .  ^ 

3,a  ~  ,aq  '  T  €a3T  — 


assume  that  q  and  t  are  linear  functions  of  6  and  e 

A 


a3 


H  Aa.  _nap 
a  =  A  0  +  B  e  _ 

,a  a3 


■^v  =  cMVae  + 


,a 


a-3  * 


Let  us  assume  that  B  and  e  n 

,a  a3 


i^U  =  t^VU  =  0  we  have 


may  be  set  arbitrarily, 


HOy  =  b^u  =  c^VaU  =  D^VapU  =  0  . 

|!  |i  p,  (i 


coordinate  system  wherein  =  (0,0, 0,1).  Then 


*  * 
ea4  "  A  =  0 


;  =  t 


*  *  0  . 


Since  e  is  symmetric  we  may  suppose  that 

a3 


gUa-3  _  Bu3a 


Dnva3  =  D^v3a 


A1.7 
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We  may  without  loss  of  generality  suppose  that 


A1.8 

=  0 

A1.9 

Qlivalt  *  ^vlKt  =  0  _ 

Let  us  now  assume  that  the  fluid  is  isotropic.  That  is,  suppose  that 

A|ia  „M.va  J-ivap  .  ,  .  _  .  .  - 

A  ,  B  ,  C  ,  D  are  invariant  under  rotations  of  the  spatial 


coordinates.  It 

is  well  known  (Jeffrey:  Cartesian  Tensors)  that  this 

requires 


A1.10 

AiJ  t  4g1J 

Bijk  5  Beljk 

Ccljk 

Dijkx  .  Digijgki  +  ^jj  +  „  ligjk 

Then  equations  A1.6  and  the  skew-symmetry  of  imply 

* 

B  =  0  . 

Since  A1.4,  A1.8  imply 


Al.ll 

BaP*  :  o 

when  a,  P,  or 

y  is  b  f  we  have 

A1.12 

Bap7  ;  o  . 

p,v 

The  equation  A1.3  and  the  symmetry  of  t 


show  that 


YT 
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A1.13 


c  t  0  . 


Equation  A1.5  shows  that 


Al.lU 


*  c4na  *  Q 


^  Q  myj  ^  ^ 

The  only  non-zero  components  of  C  '  are  those  of  the  form  C  ~ 
satisfy  isotropy  we  must  have 


cljU  2  c  gl* 


Thus 


A1.15 


CapV  =  C,  Aapuv  . 


4<x  * 

Since  A  =  0  we  have 


A1.16 


14 

1  0  0  A  i 


0  1  0  A' 
0  0  1  A 
0  0  0  0 


24 


34 


If  we  define  a  vector  A  by 


A1.17 


we  have  the  result 


=  A  ^Ve  „  + 
) v 


To 


A1.18 
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Thus  a  change  in  0  with  time  causes  a  heat  flow  in  the  direction 


M 


of  A  .  This  contradicts  isotropy  unless 


A1.19 


=  0 


We  may  now  write 


A1.20 


A^V  =  Af 


In  view  of  equations  A1  9^  AI.4,  A1.10,  we  may  put 


A1.21 


^  =  D  A^V  AaP  +  D0  A^a  AVP  +  D,  iP  AV<X 
12  5 


The  forms  of  qa  and  TaP  are  the  11 


A1.22 


q^1  =  A 


A1.23 


r^V  =  C  A^V0  +  D,  A^V  A°P?  +  (D  +D  )  A^aAVPc 

1  aP  2  y  ' 


'aP 


Observe  that  A1 . 1  implies  that 


A1 ,2k 


a 


,a 


0(c2) 


Since 


A1.25 


CL  Cl  1  (,CL  -fcOL 

s  =  spu  +  ^  £  -  ei 


We  have  in  view  of  2.5.8 
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A1.26 


Since 


A1.27 


we  have 


S  pUa  +  S(pUa)  +  0(e2)  . 

)  CL  jCL 


(nua  +  if)  „  =  0 

,a 


Ai.28 


(pUa)  =  0(e2) 

>  a 


Thus 


A1.29 


0(€2)  . 


Now 


A1.30 


be  \  ■  2, 

5  lp  +0(0  • 


According  to  A1.28 


A1.31 


b  =  -  pUa  +  0(e2)  . 

,a. 


2 

Thus,  to  within  0(e  ) 

A1 . 32  0  = 


f  80  N\  TIa  f  80 

Js'  V  ="v^ 


0  € 


a 

a 
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The  form  of  t^V 


is  then 


A1.33 


Dx  A1 


|1V 


) 


a 

e 

a 


(D2  +  D?)  AV|i 


"aP 


Li 

Substitution  of  A1.22,  A1.33>  into  Al.l  finally  gives  us  the  form  of  q 
and  t  as 


=  -K  A^V0 


K  >  0 


T^V  .  -X,  A^  Apv  e  .  -  V  A^V  Aap  e  B  ; 
1  ap  2  ap 


\  >  0  . 
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A2.1 

Thus 

A2.2 


A2. 3 

Now  P 


APPENDIX  2 

REPLACEMENT  OF  0(e2p)  +  0(e2|j)  by  0(c2P)  . 


We  may  write  the  equation  of  state  in  the  form 

p  =  p(u,p)  ;  U  =  (n-p)c2  . 


P*-P 

P 


b?  \  p*-p 

<Wu  P 


c2 
P  C 


P 


'  2  2n 
€  PC 


Now 


[1C 


2.  2  2 

>  >  e  ,  so  we  may  replace  0(€  pC— )  by  0(e  )  . 


Thus 


P*-P 

P 


0(c2) 


[1- 

i-C2 


cf.  Taub  [4]  so  we  may  write 


1^1 


ri_ 

Lc2 


A2.4 
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Under  classical  conditions  (^)  and  ~  *  (^7)  )  are  0 ( 1 ) 

^  L  c  ^  -jl 


Under  relativistic  conditions  all  substances  exist  as  gases.  It  has  been 
shown  (cf.  Israel,  Bas1'  topics  of  the  relativistic  theory  of  shock  waves) 
that  for  Eoltzmann,  Feimi,  and  Bose  gases 


A2.5 


We  may  conclude  that  for  all  "reasonable"  equations  of  state 


P 


=  0(€2) 


under  all  conditions. 
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Let 


APPENDIX  3 


PROOF  THAT  pG  =  ')  p  9  . 


Consider  a  region  of  volume  ,  containing  fluid  at  equilibrium. 


U.  =  Total  energy  in  the  region 


S.  =  Total  entropy  in  the  region 


A3.1 


M  =  Total  mass  of  fluid  A  in  the  region 


Define 


T  =  Temperature  of  the  region 


P  =  Pressure  of  the  region 


G1  =  ui  +  P1V1  -  T1S1 


It  is  well  known  (cf.  Zemansk}^  Heat  and  Thermodynamics;  McGraw  Hill, 
pp.  417-419)  that 


A3. 2 


dG  =  SdT  +  VdP  + 


N 


k=l 


^kdMk 


and 


N 

:i  ■  X  “A 

k=l 


A3. 3 
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Now  our  6 


is  related  to 


by 


A3. 4 

Q  —  ^ 

pVjT  * 

Hence 

N 

K?) 

A5.5 

0  =  /_,  (  T  , 

A=  1 

Differentiating 

A3 .4  yields,  with 

A3. 2 

A3. 5 


A3.6 


Since 


A3. 7 


We  have 


A3, 8 


If  we  then  put 


A3- 9 


ae 


OT 


dP  +  r)d 


+ 


dS  -  d  (Hr)  - de 

TdS  =  d(*)+Pd  (£)-£  »k  d 


Equations  3*2.1  and  3*2.5  follow. 
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